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1.  Introduction 


An  antenna  array  is  composed  of  many  elements  whose  excitation  am¬ 
plitude  and  phase  can  be  individually  adjusted  to  yield  a  desired  pattern. 
If  all  array  elements  operate  properly,  then  well-known  analytic  techniques 
can  be  used  to  find  the  optimum  amplitude  and  phase  of  each  element,  to 
yield  a  given  beamwidth  and  sidelobe  level.  However,  if  any  elements  fail, 
then  no  analytic  means  exists  to  find  the  aperture  distribution  that  compen¬ 
sates  for  the  degradation  of  the  pattern.  Element  failures  destroy  symmetry 
and  cause  sharp  variations  in  the  field  intensity  across  the  array  aperture, 
increasing  the  sidelobe  level  of  the  power  pattern.  The  simplest  solution  to 
this  problem  is  to  increase  the  taper  of  the  array  distribution  in  order  to 
lower  the  sidclobc.s  back  to  the  design  level.  Unfortunately,  this  solution 
has  several  disadvantages.  First,  as  the  design  sidelobe  level  for  an  array  de¬ 
creases,  the  sensitivity  of  the  sidelobe  level  to  a  perturbation  in  an  excitation 
value  increases.  Second,  changing  the  design  taper  may  broaden  the  main 
beam  to  an  unacceptable  level.  Finally,  merely  increasing  the  taper  does  not 
compensate  for  lack  of  symmetry  or  smooth  out  sharp  field  variations.  Com¬ 
pensation  for  the  defective  elements  can  be  achieved,  by  numerically  finding 
the  excitation  of  each  nondefective  element  that  optimizes  some  function. 
This  function  must  take  into  account  the  location  of  the  failed  elements. 

The  approach  used  in  this  study  is  to  reformulate  the  optimization  prob¬ 
lem  to  take  into  account  the  element  failures.  An  algorithm  is  developed  that 
synthesizes  the  amplitude  and  phase  of  each  nondefective  element  in  order  to 
produce  the  lowest  sidelobe  level  for  a  given  beamwidth.  This  algorithm  can 
also  be  used  to  find  the  clement  failure  limits  of  any  array.  More  specifically, 
the  algoriUim  can  find  the  maxiiiuim  number  of  elements  beyond  which  it 
is  impussiljle  to  recover  lost  performance.  Also,  the  loss  in  performance  can 
be  dctcrmiiu'd  by  the  geometric  arrangement  of  the  failed  elements.  The 
miiiierical  iinplemcntation  of  the  synthesis  algorithm  involves  the  minimiza- 
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tion  of  a  nonlinear  function,  which  is  the  ratio  of  the  sum  of  a  set  of  peak 
sidelobe  power  points  to  the  power  in  the  main  beam.  This  minimization 
will  be  carried  out  via  a  conjugate  gradient  method. 


II.  Power  Pattern  With  Element  Failures 


This  analysis  neglects  the  element  factor  and  considers  only  the  array 
factor.  No  mutual  couplinp  effects  are  modeled.  Failed  elements  are  deleted 
from  the  array  factor.  T\>  ”  time  convention  is  understood,  where  u>  is 

the  angular  freque^'  .  .  u  '  •  xcitation  with  a  freespace  wavelength  of  Ap. 
Letting  N  denote  ti-'  ;  mh  r  of  nonfailed  elements,  the  array  factor  (AF)  is 
given  by 

A 

AF  =  +  (1) 


kt  =  tos(0)  (2) 

ky  =  koS\t\{d)s\i\(<p)  (3) 

and  ko  =  2ff/A„.  The  excitation  of  the  nth  element,  denoted  by  i„,  can  be 
exj-'rossed  })y  the  complex  quantity  In  =  Un  + jin,  where  v'Un  •*'  yields  the 
amplitude  and  (n)i"‘(vn/n„)  yields  the  phase  angle.  The  power  pattern  is 
proi)ortii>aal  to  the  square  of  the  magnitude  of  AF. 

The  array  is  assumed  to  operate  in  either  sum  or  difference  mode.  The 
sum  mode  produces  a  peak  in  the  pattern  in  the  broadside  direction  and 
the  difference  mode  produces  a  ludl  in  the  broadside  direction.  A  possible 
clement  configuration  to  generate  both  types  of  patterns  is  shown  in  Fig.  1. 
Sum  mode  oi)eratiou  hIUavs  all  37  elements  to  be  independently  controlled, 
as  shown  in  Fig.  1(«  ).  The  same  array  operating  in  difference  nntde  has  only 
17  independently  controlled  elements,  as  shown  in  Fig.  1(b).  This  is  because 
half  of  the  array  elements  are  excited  with  their  phase  opposite  from  that  of 
(lie  lorn  sponding  <-!cinents  symmetric  with  res])ect  to  tlie  center  of  the  array. 
If  an  clement  fails,  then  the  symmetric  counterpart  is  deactivated.  This 
ensures  n  stationary  beam  null  and  a  symmetric  beam  slope,  independent  of 
the  synthesis  techniejue  u.sed. 
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(a) 


O  O  •  • 

o  o  •  • 

o  o  o  •  •  • 

o  o  o  •  •  • 

o  o  o  •  •  • 

o  o  •  # 

o  o  •  • 

(b) 


Fig.  1.  Independently  Controlled  Elements  for  (a)  Sum  Pattern  and  (b)  Differ¬ 
ence  Pattern 
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The  synthesis  of  the  excitation  values  has  three  possible  combinations: 
synthesizing  with  amplitude  only,  synthesi?'  .5  with  phase  only,  and  synthe¬ 
sizing  with  both  amplitude  and  phase.  It  is  well  known  that  synthesizing 
the  lowest  possible  sidelobe  level  for  a  given  beamwidth  over  a  symmetric 
aperture  yields  a  uniform  phase  and  a  nonuniform  amplitude  distribution. 
However,  if  the  failed  elements  yield  a  nonsymmetric  array  distribution,  then 
the  optimal  phase  may  not  be  constant.  Phase  synthesis  with  a  constant  am¬ 
plitude  was  found  by  Deford  and  Gandhi  (Ref.  1)  to  require  an  extremely 
large  iiii;nber  (if  elements  to  yield  low  sidelobes.  A  similar  analysis  by  Peters 
(Ref.  2)  found  that  phase-only  synthesis  was  not  very  effective  in  reducing 
sidelobe.':  for  flxcd-amphtudc  arrays.  Therefore,  both  amplitude  and  phase 
synthesis  would  seem  to  be  required  for  the  generation  of  a  sum  pattern. 
Since  the  difference  pattern  preserves  the  symmetry,  amplitude-only  synthe¬ 
sis  is  required.  The  sum  power  pattern  is  given  by 

,v  1  *  f  ^ 

VjunCOsfVn)  -  VnSin(v,.l]  -t  ^|u„  sin(  Vn  )  -l'  v„  cos(  Vn  )1 

’1-1  .n  r  1 

where  v,,  -  A:,?,,  +  and  the  difference  power  pattern  may  be  written  as 

'A'  1  ^ 

P{e,<p,u)=  ^u„sin(v>n)  (5) 

►nrl 

Note  tluit  the  amplitude  is  allowed  to  be  negative  on  the  difference  pattern, 
since  that  invfilves  only  a  simple  phase  shift. 


P{B,  0.  u,  ti)  =. 
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III.  Minimizing  Nonlinear  Functions 


Equations  (4)  and  (5)  indicate  that  the  power  pattern  of  the  array  is  a 
nonlinear  function  of  the  amplitude  anti  phase  of  the  excitation.  Minimizing 
the  peak  sidclobe  level  will  require  a  mathematical  technique  that  can  handle 
an  nonlinear  function.  A  multivariable  quadratic  function  can  be  minimized 
by  the  conjugate  gradient  algorithm  developed  by  Hestenes  and  Stiefel  (Ref. 
3).  An  arbitrary  nonlinear  function  can  also  be  jninimize^by  the  same  algo¬ 
rithm  if  a  suitable  quadratic  approximation  can  be  found.  Although  several 
variation.s  are  given  in  the  literature,  all  have  the  same  two  major  steps  d:s- 
c\issed  by  Hestenes  (Ref.  4).  First,  the  function  must  be  approximated  by  a 
quadratic  truneation  of  a  Taylor  scries  expanded  around  an  initial  estimate 
of  the  sohition.  Second,  this  quadratic  is  minimized  by  the  conjugate  gradi¬ 
ent  algorithm  to  <d)tain  a  new  estimate.  This  new  estimate  is  in  turn  used  as 
the  expansion  point  for  a  new  quadratic  approximation.  Since  the  quadratic 
approximation  is  valid  only  near  the  point  of  expansion,  finding  a  minimum 
of  this  (luadratic  is  the  same  as  finding  a  minimum  of  the  original  function 
only  if  botli  have  the  same  minimum.  Therefore,  the  algorithm  must  be 
restarted  so  that  the  quadratic  approximation  can  ultimately  be  expanded 
around  a  point  that  approaches  the  minimum  of  the  original  function.  The 
one-diiiunsional  ca.se'  is  illustrated  m  Fig.  2.  Given  the  function  /(«)  and  the 
initial  gtiess  the  quadratic  qo(«)  is  formed  and  given  by 


Ocl.i)  -  fiSu)  ■*  /*(»(.)(»  «o) 


^  f"\ «f>)(  «  - 
2 . 


(6) 


The  mininnmi  of  this  quadratic  is  denoted  by  >i  and  is  given  by  the  Newton 
estimate 


Tile  estimate  «,  is  then  ii.scd  t<»  com])uto  <7,(0).  The  iteration  then  contin¬ 
ues  until  the  niininmiii  of  the  quadratic  nppro.xiinntion  coincides  with  the 
iiiinimuiu  of  the  original  function. 
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Fig.  2.  One- Dimensional  Minimization 
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The  principles  of  the  one- dimensional  minimization  apply  also  to  mul¬ 
tidimensional  functions.  For  example,  consider  the  function  /(u,v),  where  u 
and  f  are  N  x  i  column  vector  variables.  The  quadratic  approximation  may 
be  written  in  matrix  form  as 

f(u,v]  _  +  h-7'rG+  (8) 

where  T  denotes  the  transpose  and  the  components  are  defined  as  follows. 
The  vector  7  is  a  2A’  x  1  column  vector  defined  by  7  =  [u  t]^.  The  gradient, 
denoted  by  G',  is  a  '2N  x  1  column  vector  evaluated  at  and  is  given  by 

C  G\]''^,  where  each  subvector  has  the  form 


(9) 


The  Hessian,  denoted  by  H ,  is  a  2N  x2A'  matrix  whose  elements  are  evaluated 
at  (v'.r')  and  exj)ressed  as 


H 


Huu 


(10) 


where  cacli  submatrix  has  the  form 


r  -Plf. 

dxiidv  1 

-  B’J 
flu |flV} 

. .  -J— 

duyOv  N 

=■ 

gy 

flu  jflr  I 

flvjflrj 

..tlL 

flu  >fll  n 

a’; 

flu  fvflv  1 

flu  ^flv  ) 

dundv  s 

I'sing  the  one-diniensional  analogy,  the  minimum  of  Eqn.  (8)  may  be  coin- 
l)uted  from  the  matrix  form  of  Newton’s  formula,  as  given  by 


7i  =  7o  + //  ‘(7o)G’(7o)  (12) 

However,  this  recpiires  that  the  inverse  of  H  be  c(.>mputed.  In  many  cases 
of  jjractieal  intiuest.  only  an  estimate  of  the  minimum  of  the  quadratic  is 
retjuired;  therefore,  it  is  more  efficient  to  use  the  conjugate  gradient  method 
VO  obtain  .'in  estimate  (.if  this  minimum.  Making  the  substitution  =  u  ~  u\ 
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2,  =  t;  -  i’,  and  D  =  ~G  into  Eqn.  (8)  and  excluding  the  constant  yields  a 
standard  quadratic  function,  denoted  by  c  and  written  as 


Minimizing  this  expression  is  the  same  as  solving  the  linear  system 


(13) 


(14) 


where  the  Hessian  matrix  H  is  always  symmetric  and  nonsingular  but  not 
always  positive  definite.  Therefore,  the  conjugate  gradient  algorithm  must 
be  general  enough  to  handle  this  situation.  A  suitable  algorithm  w'ith  L 
restarts  and  h  <<  A’  conjugate  gradient  iterations  per  restart  is  given  as 
follows  (Ref  5); 


begin  loop  for  h: L 

compute  D' ,  //'  and  tel  =0,2^  =0 

r'  =  D' 

U.V  *^U,9 
(Jxi  — 

i - - \ 

(7«.9u/  -»  Wv.^W 

Pu,v  ~  0OQu,t 

begin  loop  for 

-  //uupi  + 

—  ^IvuPu 

1 

'»*  -  ^ - — : - ; 

^i,V‘  iJl.r  +  ^kPu.v 

^u.r  ~  ^u.v  ‘**9u.r 

t’-J.’-i)  +  (’■J.’-J) 


(16) 

(16) 

(H) 

(18) 

(19) 


(20) 

(21) 

(22) 

(23) 

(21) 

(26) 
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if  err  <  tolerance,  terminate  k  loop,  else 


-continue  k  loop 


(26) 

(27) 

^  1 

(qu.qu)  +  (9.,9.> 

(28) 

pi.V  =  pi,.  +  0kqu,t 

(29) 

U,v'  +  '  =  U.v'  -t-  2^  +' 

(30) 

continue  I  loop 

The  I  loop  determines  the  quadratic  approximation  to  the  function  f{u,v). 
The  gradient  and  Hessian  are  computed  at  the  current  estimate  and  the 
ai)propriate  conjugate  gradient  parameters  are  initialized.  The  k  loop  obtains 
an  estimate  of  the  minimum  of  c  by  computing  an  estimate  of  the  solution 
to  the  linear  system  described  by  Eqn.  (14).  The  dominant  computational 
effort  for  each  k  iteration  is  caused  by  the  matrix  vector  products,  which 
require  CP(8;V*)  imiltiplications.  The  computation  of  the  Hessian  requires 
C>(tTiV^)  multiplications,  where  T  is  the  number  of  multiplications  per  matrix 
element.  Since  T  depends  on  the  function,  the  computational  efficiency  of 
the  overall  algorithm  is  function-dependent.  An  alternate  algorithm  used  by 
Fong  and  Birgenheier  (Ref.  6)  and  Press  et  al.  (Ref.  7)  does  not  require 
the  ccunputation  of  the  Hessian,  but  implicitly  assumes  that  it  is  positive 
definite.  This  can  lead  to  an  unpredictable  breakdown  in  the  algorithm. 
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IV.  Sum  Pattern  Synthesis 


The  goal  of  the  sum  pattern  synthesis  is  to  maximize  the  power  in  the 
main  beam  direction  while  simultaneously  minimizing  the  power  in  the  peak 
sidelobe.  Unfortunately,  forcing  the  sidelobe  to  a  lower  level  et  a  particular 
angular  location  without  any  constraint  on  the  other  sidelobes  will  merely 
shift  the  peak  sidelobe  to  a  new  location.  Because  the  peak  sidelobe  location 
may  change,  there  is  no  mathematically  elegant  way  to  achieve  the  goal  di- 
rect'y.  Rather,  the  same  objective  may  be  achieved  indirectly,  by  nrinimizing 
the  average  of  the  accumulated  peak  sidelobe  power  points.  The  average  is 
defined  to  he  the  average  of  these  peak  sidelobes,  computed  at  each  restart 
of  the  minimization  algorithm.  If  Mt  denotes  the  number  of  distinct  peak 
sidelobe  points  at  the  /th  restart,  then  the  average  power,  denoted  by  Pa{u,v), 
is  defined  by 

1 

/’.(ti.v)  =  —  (31) 

The  jjower  at  tlie  ith  angular  point  is  defined  as  P,  =  P(u,v,0i,(p,),  where  i  =  0 
denotes  the  main  beam  location  and  i  =  1...M  denotes  peak  sidelobe  points. 
In  order  to  minimize  the  sidelobe  power  while  simultaneously  maximizing 
the  main  l)eam  power,  a  function  /(«,v)  is  defined  as  the  ratio  of  the  average 
I)eak  i)owei  to  the  main  beam  power  and  is  expressed  as 


/(v.t’) 


Pa{u,v} 
Po(“, v) 


(32) 


Minimizing  this  function  will  acliicve  the  desired  goal.  A  possible  one- 
dimcnsioiiul  scenario  is  illustrated  in  Figure  3.  The  peak  sidelobe  level  is 
found  at  tlie  initial  start.  This  sidelobe  level  is  then  reduced  by  a  partial 
minimization  <A  f[u,v).  At  the  first  restart  the  peak  sidelobe  point  is  still  at 
the  same  lorntion.  so  that  peak  is  reduced  again.  At  the  second  restart  a 
new  jK-nk  sidelobe  location  is  found,  so  this  j)cak  and  the  previous  peak  are 
averaged  and  reduced.  At  the  third  restart  tlie  peak  location  has  changed 
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again,  so  the  average  is  updated  and  /{«,«)  is  minimized.  At  the  fourth 
restart  another  point  is  added  to  the  average  and  f{u,v)  is  minimized.  The 
pattern  has  reached  the  steady-state  power  distribution. 

An  inherent  characteristic  of  this  method  is  that  all  the  sidelobes  will 
tend  to  be  brought  to  the  same  level.  It  should  be  noted  that  all  the  sidelobe 
sample  points  could  have  been  used  a  priori  at  each  restart  without  searching 
for  the  peak  sidelobe  point  and  the  results  would  be  the  same.  However,  using 
the  accumulated  average  has  several  advantages.  The  function  /(u,  v)  changes 
rapidly  when  few  sidelobe  points  are  present,  so  that  the  peak  may  be  brought 
down  fast.  As  more  peak  points  are  added  to  the  average,  a  self-damping 
effect  takes  place  and  reduces  the  sensitivity  of  the  pattern  to  changes  in 
u  and  f,  allowing  a  smooth  transition  to  an  optimum.  Accumulating  the 
peak  points  is  numerically  much  faster  than  starting  with  all  points,  since 
the  computation  time  for  evaluating  /(u,  v)  is  dependent  on  A/j. 

The  sidelobe  region  is  defined  as  the  hemisphere  formed  by  the  limits 
0  <  <  2-  and  0  <  0  <  rr  minus  the  circular  aperture  of  the  main  beam  up 

to  the  first  null.  The  density  of  the  sample  points  should  be  chosen  in  such 
a  way  that  the  peak  sidelobe  may  be  accurately  estimated.  A  good  rule 
of  thumb  is  that  the  sample  separation  distance  should  be  no  greater  than 
half  the  main  beamvvidth  of  the  design  pattern.  A  uniformly  spaced  square 
lattice  is  sufficient  and  is  shown  in  Fig.  4.  The  parameter  is  the  radius 
of  the  main  beam  null  of  the  design  pattern,  increased  slightly  to  allow  the 
beam  to  broaden  as  the  sidelobe  level  is  reduced. 

The  components  of  D  and  //  required  for  the  conjugate  gradient  method 
are  given  as  follows: 


8f 

1 

Pa  dPn 

8  Pa 

Po 

.  Po  8um 

5u„ 

/?;;*  -- 

8f 

I 

Pa  8Po 

8l,n 

Po 

.  Po  8v^ 

8v,n  \ 

II"'’'  .  — - 

uu  a 

OU 


du„ 


1 

Po 


nPp 

lifTl 


SFo 


8^  Pa 


8u„  du„8u„ 


Pa  d^Po 
Pn  dum8u„ 


(33) 

(34) 


(33) 
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«mn  ^  gy  _  ±  [nn  ^  ^  g^^O 

Po  r*  «u„  “  5v„  ^  8u^dv„  Pogu„«v„ 


=  ^^r,"  (37) 

gy  Mzj"^^°  I  I  ^38) 

"  Pol  ’gtrm  'dv„^ev„dv„  Podv„dv„\  ^  ^ 

The  power  at  the  «th  angle  can  be  computed  by  introducing  the  auxiliary 
functions  and  Xi  S’Hd  defining  Pi  =  +  Xi  such  that 


=  5Z(**“  coafV-tti)  --v„  sinlVitti)] 


Xi  -  sin(t/'«,)  +  v„  coslV-wi)]  (40) 

to  =1 

where  vv,  -  +  fiyiVu-  The  first -derivatives  of  the  power  used  in  the 

gradient  are  given  by 

d  P 

— *-  =  2[fiCOs(vmi)  +  Xii>ia(Vmi))  (41) 

dP 

=  2(x<cos(v„<) -(isinfVmi)]  (42) 

and  the  second  derivative  terms  used  in  the  Hessian  are  expressed  as 


O  -a  -  =  Scosltymi  "  V’ni) 

du^dUn 
8^  Pi 

- - -  =  2sin(V»mi  -  V'ni) 

8^  Pi  _  8^  Pi 

8v„dv„  5u„5u„ 


Note  that  the  computation  of  the  Hessian  requires  the  evaluation  of  the 
power  at  the  main  beam  location  and  at  each  peak  sidelobe  point.  Since 
the  power  computation  requires  0(N)  multiplications,  the  multiplications 
per  matrix  element  is  T  =  O(MiN).  Therefore,  the  Hessian  computation 
dominates  the  computation  of  a  single  I  iteration. 
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V.  Difference  Pattern  Synthesis 


The  synthesis  procedure  for  the  difference  pattern  is  similar  to  that  used 
for  the  sum  pattern,  except  that  now  only  amplitude  synthesis  is  necessary. 
The  appropriate  function  defined  in  Eqn.  (32)  is  reduced  to  one  dimension 
as 


/(tt) 


Po(u)- 


The  jiower  al  the  ith  location  is  defined  as  P  =  where 


N 

^  u„  sin(vu.i) 

«=1 

The  first  derivative  of  the  power  is  computed  from 


(46) 


(47) 


— -  2£<  sin(V>„. ) 

V  tifn 


with  the  second  derivative  given  by 


d'Pj 

dULmdUn 


2sin(vmi)sin(Vni). 


(48) 


(49) 


Fvi  the  configuration  shown  in  Fig.  1,  the  power  pattern  is  symmetric  across 
the  y  -  z  plane.  Thus,  the  fc-space  sample  region  need  only  contain  points 
where  k,  >  0.  as  shown  in  Fig.  5. 


VI.  Results  and  Discussion 


A  possible  element  failure  scenario  is  shown  in  Fig.  6.  The  hexagonal 
array  has  a  total  of  91  elements  spaced  O.6A0  apart  with  3  element  failures. 
The  desired  aperture  function  for  the  sum  pattern  is  a  parabolic  with  power 
N  -  2  and  an  edge  taper  of  -20  dB.  This  yields  a  -35  db  sidelobe  level.  The 
aperture  function  for  the  difference  pattern  is  a  Bayliss  distribution  with 
ft  =10  and  a  -40  dB  sidelobe  level.  All  the  power  patterns  are  plotted  in  dB 
using  cylindrical  coordinates,  with  $  as  the  radial  coordinate  in  the  range 
0<  6  <90=.  A  reference  -32  dB  floor  is  placed  on  each  plot. 

The  siiiTi  j)owcr  pattern  with  no  element  failures  is  shown  in  Fig.  7. 
The  peak  sidelobe  level  is  -30.9  dB  with  a  peak  gain  of  20.8  dB.  The  main 
beam  has  azimuthal  symmetry  with  a  half-power  beamwidth  of  12.86°.  The 
synthesized  sum  pattern  of  Fig.  8  required  300  restarts,  with  5  conjugate 
gradient  iterations  per  restart,  to  solve  for  the  176  required  unknowns.  The 
effect  of  the  failures  is  to  increase  the  sidelobe  level  to  -27.4  dB  and  reduce  the 
gain  to  20.5  dB.  The  main  beam  becomes  slightly  elliptical  with  a  minimum 
half-power  beamwidth  of  12.57°  and  a  maximum  of  13.00°.  The  synthesized 
pattern,  shown  in  Fig  9,  has  a  peak  sidelobe  level  of  -31.3  dB  and  a  gain  of 
19.5  dB.  The  main  beam  also  has  an  elliptical  cross  section  with  a  minimum 
half-power  beamwidth  of  13.71°  and  a  maximum  of  16.29°.  As  expected, 
the  reduction  in  the  sidelobe  level  produces  a  corresponding  increase  in  the 
beamwidth  and  a  drop  in  gain. 

The  difference  power  pattern  with  no  element  failures  is  shown  in  Fig. 
10.  The  i>eak  sidelobe  level  is  -33.0  dB  with  a  peak  gain  of  18.4  dB.  The 
angular  half-power  null  width  between  the  peaks  is  8.57°.  Fig.  11  shows  that 
(he  effect  f»f  the  element  failures  is  to  increase  the  sidelobe  level  to  -21.7  dB 
and  dro])  the  gain  to  18.0  dB,  with  the  half-power  null  width  remaining  about 
the  same.  The  synthesized  difference  pattern  of  Fig.  12  required  100  restarts 
with  5  conjugate  gradient  iterations  per  restart  to  solve  for  the  80  unknowns. 
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Fig.  G.  Hexagonal  91  Element  Array  with  3  Element  Failures 


Avi  th  N  u  Failu  res 


Fig.  i).  Syntlu'sizcd  Sum  Puwfr  Pattern  with  3  Element  Failures 
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Fig.  10,  IDiffereiice  Pi/wci  Pattern  with  Nu  Failure 


The  peak  sidelobe  level  is  -28.7  dB  with  a  gain  of  17.1  dB.  The  sidelobe  level 
of  the  synthesized  pattern  cannot  be  brought  down  to  the  level  of  the  design 
pattern  as  was  the  case  with  the  sum  pattern,  because  the  degrees  of  freedom 
are  now  about  halved.  The  half-power  null  width  between  peaks  increased 
to  9.09'. 

The  results  of  this  test  case  and  several  others,  indicate  that  it  is  possible 
to  reduce  the  sidelobe  level,  of  an  array  with  failed  elements,  down  to  near 
the  design  level.  However,  the  success  depends  on  the  nature  of  the  failure. 
For  example,  if  the  central  element  fails  in  a  sum  pattern,  the  synthesis  yields 
no  improvement.  This  is  to  be  expected  since  this  failure  is  quite  similar  to 
blockage.  As  expected,  elements  which  fail  and  have  a  large  relative  weight 
are  much  more  difficult  to  recover  from  than  elements  with  a  small  weight. 
Also,  clustered  failures  are  easier  to  recover  from  than  the  same  number  of 
random  failures.  This  is  because  a  clustered  failure  yields  a  relatively  local 
increase  in  the  sidelobe  region  as  opposed  to  the  same  number  of  random 
failures  which  yields  a  uniform  increase  in  the  sidelobe  level.  It  was  iriterest- 
ing  to  observe  that  is  was  easier  to  recover  from  an  entire  row  failure  than 
from  the  same  number  of  random  failures.  In  general,  the  decrease  in  the 
sidelobe  level  comes  w'ith  a  drop  in  gain  and  a  slightly  broader  main  beam. 
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VII.  Conclusion 


The  performance  degradation  of  arrays  with  element  failures  may  be 
partially  compensated  for  by  a  redistribution  of  the  amplitude  and  phase  over 
the  remaining  elements.  The  extent  of  this  compensation  is  determined  by 
the  number  and  location  of  the  failed  elements.  The  accumulated  averaging 
scheme,  combined  w’ith  the  conjugate  gradient  algorithm,  provides  a  very 
_stable_means  to  synthesize  the.new.  array  .aperture,  distribution.  _ 
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Appendix:  Fortran  77  Programs 


This  appendix  contains  the  Fortran  77  programs  SOHSTN  and  DZFSTN  as  well 
as  all  associated  subroutines.  The  program  SOMSTN  performs  the  sum  pattern 
synthesis  with  element  failures  and  the  program  DIFSTM  performs  difference 
pattern  synthesis  with  element  failures.  The  required  inputs  to  each  program 
are  listed  in  the  comments  found  in  the  source  code.  These  programs  were 
used  to  generate  all  the  results  presented  in  this  report. 
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1 

1  c 

3  c 

4  C 
6  C 

6  C 

7  C 

8  C 
e  C 

10  C 
n  C 
13  c  - 

13  C 

H  c 

16  C 
18  C 

17  C 

18  C 
18  C 

30  C 

31  C 
22  C 
33  C 
24  C 
35  C 
2fl  C 

27  C 

28  C 
28  C 

30  C 

31  C 
33  C 

33  C 

34  C 

35  C 
30  C 

37  C 

38  PARAMETER  (NL=6 ,NTE=91 , NBE=3 ,NE=e8 ,NSP=1660 ) 

39  REAL*4  0(NE).V(NE),I(IIE).T(KE),SU(RE),SV(»E),CU(»E),GV(NE) 

40  REAL*4  HUO(HE,HE),H0V(HE,HE).H7U(IE,irE),HVV(HE,HE) 

41  REAL*4  PU(NE)  ,PV(1TE)  .QO(RE)  .QV(RE)  ,RTI(HE)  .RT(HE)  .2U{RE)  .ZV(HE) 

42  REAL*4  XX(NSP)  ,KT(IISP) 

43  IHTECER  IP(NSP) ,BAD(NE) 

44  0PEN(DNIT=2,FILE= ’OUTPUT’ ,STATOS=’0HKHOWN’) 

45  0PEN(nNIT=4,FILE=’S0M_FHED_EST’ ,STATUS= 'UNKNOWN ’ ) 

46  C  ********♦*•»*•••****•*♦•*«••••*•••*****»•***♦*•*♦•••***•**••*♦***** 

47  C  *  GENERATE  THE  GEOMETRT .  * 

48  C  •*»***••*•••»•••********»«•**••**♦****♦***•♦•*•»••**♦**♦*****•**•** 

49  DS=0.6 


PROGRAM  SUHSTN 


THIS  PROGRAM  STHTHESIZES  THE  AMPLITUDE  AND  PHASE  DISTRIBUTION 
NECESSARY  TO  GENERATE  A  SUM  PATTERN  OVER  A  HEXAGONAL  ARRAY  OF 
POINT  SOURCES  WITH  SOME  ELEMENT  FAILURES. 


TIMOTHY  J.  PETERS  LAST  UPDATED 

THE  AEROSPACE  CORPORATION  3/1/01 

2360  EAST  EL  SEGUNDO  BOULEVARD. 

EL  SEGUNDO,  CA  00246 


-INPUTS: -  -  -  -  -  -  .  -  _  .  - 

NL  -  NUMBER  OF  LAYERS  FORMING  THE  HEXAGONAL  LATTICE 
NTE  -  NUMBER  OF  TOTAL  ELEMENTS  =  l+3*NL*(NL4l) 

NBE  -  NUMBER  OF  BAD  ELEMENTS 

BAD(1..NBE)  -  INTEGER  ARRAY  HOLDING  THE  NUMBER  POSITION  OF  EACH 
BAD  ELEMENT  IN  THE  ARRAY. 

DS  -  TEE  ELEMENT  SPACING  IN  WAVELENGTHS. 

NE  -  NUMBER  OF  ACTUAL  ELEMENTS  =  NTE-NBE. 

ET  -  THE  EDGE  TAPER  OF  THE  PARABOLIC  DISTRIBUTION. 

N  -  THE  ORDER  OF  THE  PARABOLIC  DISTRIBUTION. 

BNF  -  THE  THEORETICAL  BEAM  WIDTH  NULL  FACTOR  FOR  THE  CHOSEN 
DISTRIBUTION. 

BBF  -  THE  BEAM  WIDTH  BROADENING  FACTOR.  THIS  ALLOWS  THE  BEAM 

TO  BROADEN  BEYOND  THE  DESIGN  VALUE  IN  ORDER  TO  REDUCE  THE 
SIDELOBES  WHEN  ELEMENTS  FAIL. 

NSP  -  NUMBER  OF  SAMPLE  POINTS  IN  THE  SIDELOBE  REGION. 

NS  -  NUMBER  OF  RESTARTS. 

HI  -  NUMBER  OF  CONJUGATE  GRADIENT  ITERATIONS  PER  RESTART. 
OUTPUTS . 

U(NE)  -  REAL  PART  OF  THE  EXCITATION  OF  EACH  ELEMENT. 

V(NE)  -  IMAGINARY  PART  OF  TEE  EXCITATION  OF  EACH  ELEMENT. 
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80  BAD(1)=3E 

51  BAD(2)=37 

62  BiD(3)=44 

53  CALL  GEONEKNL.DS.NB.NBE.BAD.I.T) 

54  C  **************•****•♦**♦*♦•*♦♦•♦•**••♦*•*♦*•♦••♦*•***•**••**♦**••** 

56  C  *  INITIALIZE  THE  TAPER  DISTRIBOTIOH .  * 

69  C 

67  ET=0.1 

59  N^2 

69  BNP=1.B 

90  IFLAG1=0 

91  IF  (IFLAGl  .EQ.-l)  THEN 

82  CALL  TAPER(HE,n,V,I.7.ET,H,HL,DS) 

93  REWIND  4 

94  WRITE(4,101)  (I,n(l).V(l),l=l,BE) 

85  ELSE 

09  END  IF 

97  REWIND  4 

98  READ(4.103)  (I,U(I),V<I),I=1.NE) 

96  C 

70  C  *  COMPOTE  THE  GAIN  (DIRECTIVITY)  OF  THE  ARRAT .  * 

72  CALL  GAIN(NE,n.V,I,T,CDB) 

73  C  ******••**•**♦***♦*♦***♦»•••*••  •••*******************»******ili»*»**» 

74  C  *  GET  THE  BEAM  WIDTH  BETWEEN  FIRST  NOLLS.  * 

75  C  . . . . . . . 

79  CALL  BEAH(NE,0,V,X,T,BWrN) 

77  C  *.♦*******.*****♦»**•***.*•*•••*»•♦••**•••****•*♦*••*•*••*****••••* 

78  C  •  GENERATE  THE  K-SPACE  SAMPLE  POINTS.  • 

70  C 

80  BBF=7.0 

81  CALL  KSPACE(NMP,NL,DS,BNF,BBF,KI,KT) 

87  C  •«,**•***•*••***♦***•**•*•*••••••**••••*••***•**••••***•**♦*»•••*•• 

83  C  •  COMPOTE  POWER  IN  MAIN  BEAM.  * 

9,  C 

85  CALL  P0WI(NE,D,V,1,T,0.0,0.0,P0) 

r.n  C  •»»•.♦****♦«•*»*•♦*****♦***•♦**•*****•♦*♦**♦♦♦•♦•♦•♦••****♦*••**♦*♦ 

87  C  *  COMPOTE  PEAK  SIDELOBE  POWER.  * 

88  C  •••**•••»•♦••*•**»*♦*••***•♦•••••*•••*••«*•♦•**»**•*****»**»♦*•♦♦♦♦ 

80  NP  =  0 

90  CALL  PEAKSL(HE,0,V,X.T,NSP,HMP.KI,KT,IP,NP,PMAI) 

91  C  «•»«**.•.•»»*•*♦*.*♦•«•••**•••*****•**•••*••***•*•*•••.•*••*•***•** 

P2  C  •  COMPOTE  NORMALIZED  PEAK  SIDELOBE  POWER  IN  DB .  * 

93  C  *».*•».*••***.*••***•»**»•♦»♦•**•»••••••*••»*»*•••*•*»••*•*»*♦*•**• 

04  PMAIDB=10.0*AL0G10(PMAI/P0) 

96  C 

00  C  *  PERFORM  CONJOGATE  GRADIENT  STEPS.  * 

07  C 

98  NS=300 


43 


se 

100 

10) 

lOJ 

103 

104  C 
lOG  C 
loe  C 

107 

108 
100  C 

110  C 

111  C 
117 

1  13  C 

114  C 

115  C 
1  ;6 

117  C 

118  C 
110  C 
120  100 
!’l  101 
132 

173  C 

124 

125  C 
12«  C 

177  C 
128 
120 
l.in 

131  C 

132  C 

133  C 

134 
I  35 

136 

137 

138 

139 
1  40 

141 

142 

143 
1  44 
145 
1  18 
147  1 


NI  =  3 

CALL  CCR1D(HE,KSP,1(MP,KS,HI.O,V,I.T,SO.SV,BP,IP.KI,KT,00,GV,HOU 
»,HU7,HVTJ,HVV,PU.P7,QU,QV.Ra.RV,Z0.2T) 

REWIND  4 

WRITE(4,101)  (I. 0(1), 7(1), 1=1, BE) 


•  POWER  PATTERN  GRAPHICS  OOTPHT 


REWIND  4 

READ(4,101)  (I, 0(1). 7(1), 1=1. BE) 


*  PRINT  THE  POWER  PATTERN  18  CTLIBDRICAL  COORDINATES.  * 


CALL  P0WPAT(NE.tI,7,I.T) 


*  PRINT  THE  POWER  PATTERN  IN  SPHERICAL  COORDINATES.  * 


CALL  SPOWPAT(NE,D, 7.1.1) 


*  FORMATS .  * 


F0RMAT(I6,1I,F16.B,1X.F1B.B) 

F0RMAT(IB.1X,F16.7.1X.F1B.7) 

END 

SOBROUTINE  GEOMET( NL , S , NE .NBE .BAD . I .1) 


♦  COMPUTE  THE  POSITION  OF  EACH  ELEMENT  IN  A  HEXOOOSAL  ARRAT.  • 


REAL*4  X(NE)  .T(HE) 
INTEGER  BAD(NE) .IFLAG 
HGT=(SqRT(3 .0)/2.0)*S 


*  CENTER  ROW. 


* 


K^O 

L^O 

XMIN=-NL*S 
DO  1  I=0,2*NL 
K=K  +  1 
L=L+1 

CALL  CHECK(L.NE.NBF.. BAD. IFLAG) 
IF  (IFLAG  .EQ.  1)  THEN 
K=K-1 
ELSE 

X(K)-JMIM+I*S 
7(K)=:0.0 
END  IF 
CONTINUE 
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MB  C 

M9  C 

♦  TOP  ROWS. 

ISO  C 

161 

INnK-2*nL+l 

163 

DO  2  J=1,»L 

1S3 

ihom=ij:um-i 

1  54 

IKIN=-NL*S+J*S/2.0 

I  65 

DO  3  I=0,IHDM-1 

t&e 

K=K+1 

167 

L=L+1 

I  SB 

CALL  CBBCK(L,RE,NBE,BAD.1FLAC) 

160 

IF  (IFLAC  .EQ.  1)  THEH 

100 

K=K-1 

lei 

ELSE 

182 

X(K)=IMIN+I*S 

103 

T(K)=HGT*J 

164 

END  IF 

185  3 

CONTINUE 

lee  2 

CONTINUE 

lfi7  C 

188  C 

*  BOTTOM  ROWS. 

• 

lee  C 

170 

INUM-2*NL+1 

171 

DO  4  J=1,NL 

172 

IN0M=1NUM-1 

173 

imN=-NL*S+J*S/2.0 

174 

DO  B  I=0,INUM-1 

17S 

K=K+1 

176 

L  =  L+1 

177 

CALL  CHECK (L,NE,NBE,B AO. IFLAG) 

178 

IF  (IFLAC  .EQ.  1)  THEN 

170 

K=K-1 

IftO 

ELSE 

181 

X(K)=IMIN+I*S 

182 

T(K)=-HCT*J 

IH3 

END  IF 

184  6 

CONTINUE 

185  4 

CONTINUE 

1  88 

RETURN 

187 

END 

188  C 

180 

SUBROUTINE  CHECK(L ,NE , NBE .BAD , IFLAG) 

100 

INTEGER  BAD(NE) 

101 

IFLAG=0 

19? 

DO  1  I:=1,NBE 

193 

IF  (L  .EQ.  BAD(I))  THEN 

194 

IFLAG=1 

195 

ELSE 

;  00 

END  IF 
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107  1  COHTINDE 

108  RETURN 

101  END 

100  C 

201  SUBROUTINE  BEAM(NE ,U ,7 .X.T .BWFK) 

202  C  ♦•••♦**♦♦••*♦*♦**************♦*••••*•♦*••*•♦••**••••**•***•••♦*•**• 

203  C  ♦  COMPOTE  THE  BEAM  WIDTH  BETWEEN  FIRST  HULLS.  • 

204  C  *♦**•***♦ *******•*******•>♦*••*••••••••••••*•••*••••••* •*••*••**•*** 

205  REAL  U(NE),V(NE).XIHE).T(HE).KXI,KTI,K1P,KTP 

208  RAD=.17453203E-01 

207  TP=.628318BE+01 

208  TMIN=0.0 

209  TMAX=10.0 

210  P-0.0 

211  NT=140 

212  DT=(TMAX-TMIH)/HT 

213  C  •*••*•*••*••**•*•••*•**•••*•** •••♦•••**»**»**»4i***»****«*********** 

214  C  *  FIND  THE  MAXIMUM  VALUE.  * 

215  C  ••**•*•***»•*••**•******•*••••*••••••*•****••***•*•*•••****••**•*** 

216  CALL  POWERCNE.U.V, X.T, 0.0, 0.0, PEAK) 

217  C  »*******+*»************************#***4*****»*«i******************* 

218  C  ♦  RECOMPUTE  AND  NOW  FIND  THE  POINT  WHICH  IS  AT  THE  SLL  BELOW  THE  • 

219  C  •  PEAK.  • 

220  C 

221  DO  1  K  =  l,90 

222  P  =  P+1.0 

223  DO  2  1  =  0, NT 

224  T=TMIH+I*DT 

225  TPS=TP*SIN(RAD*T) 

226  KXI=TPS*COS(RAD*P) 

227  KTI=TPS*SIN(RAD*P) 

228  CALL  P0UI(NE,D,V,I.T.KXI,KT1,P0W) 

229  IF  (POW/PEAK  .LE.  0.6}  THEN 

230  URITE(2.100)  T,lO.O»ALOC(?OW/PEAK) 

231  GO  TO  98 

232  ELSE 

233  END  I? 

234  2  CONTINUE 

235  89  CONTINUE 

236  BWFN=2.0»T 

237  WRITE(‘,*)  ’P  BWFN  '.P.BWFN 

238  1  CONTINUE 

2.39  100  F0RMAT(F10,B.1X,F10.6) 

240  RETURN 

241  END 

242  C 

243  SUBROUTINE  KSPACE( NMP , NL ,DS ,BNF,DBF ,KI ,KT) 

244  C  **»*»»*•«*.••*•»**•*••*»*•»•••••*•**•***•***•*•**♦•»*•••**•*♦••**•» 

245  C  *  THIS  SUBROUTINE  SAMPLES  THE  KX  >0  REGION  OF  K  SPACE.  • 
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2^e  C 

247 

248 
749 

250 

251  C 

252  C 

253  C 

254 

255  C 

256  C 

257  C 
25?  C 
:i60 

200  C 
291  C 
2(12  C 
2«3  C 

364 

305  C 
206  C 
20/  C 

208  C 

209 

2  70  C 

27t  C 

272  C 

773 

'!74 

275 

2  76 

277 

t:fs 

270 

280 

201 
282 
283 
?«1 

285 

286 
287 
788 
289 
2t-  ! 

?(»} 

2112 

■:t*3  2 
2W1  1 


REAI,?4  KI(»),KT(*)  .KC^JCR.m.KTT.KXMIN.KIMAI.KTMIN.KTMAI 

RAD=  .174r.3293E-01 

PI^.3141693E+01 

TP“.02a3ieBE+01 

♦  COMPUIK  THE  APPROXIRITE  DIAMETER  OF  THE  APERTURE. 

*4i4i**«**it  *«»****»4*****««****«************«#** 

D-2.0*D5*VL 


♦  COy.POTE  THE  BEAM  U.TDTH  BETWEEN  FIRST  HULLS  OF  THE  DESICN 

•  PATTERN .  ■ 


BWF.V=EIIF*2 . 439272/E 


♦  COMPUTE  THE  SMALLEST  BEAM  WIDTH  BETWEEN  FIRST  HULLS  ALLOWED  BT 

•  THE  STMTRESIS  ALGORITHM. 


BWFNS=BBF»BWFH 


*  COMPUTE  THE  K  SPACE  RADIOS  OF  THE  ALLOWED  MAIN  BEAM. 

«*********«***4*****««4*4»4**44*««.*»4*****»*******»»4*i«**4**«*»*4* 

KC=TP*0.34 

KC=8UF.'JS/2.0 


•  GENFRATE  THE  RECTANGULAR  LATTICE  OF  SAMPLE  POINTS. 

♦  ••♦•♦•♦♦•♦•♦'•♦•**i»*«4*********»»4'*****»»»*»*******»«*****»****** 

KTMin---TP 

KT.MAI-TP 

KIMIK=-Tr 

KIMAl=Tr 

NKX^4D 

HK7^4& 

DKX-^(KXMAX -KIMIN)/(NKI-  t) 

CKY=(KTMAX-KIMIN)/(NKT-1) 

NMP-0 

HO  1  I  =  y,I!KT 

KTT-K7KIN'+1«DKT 
DO  2  J  O.NKI 
EXX^.KIMINN  I-DKX 
Kn=5QhT(KII*KrX-tKTy«KTT) 

IF  (ly.R  .GT.  KC)  .and.  (KR  .L7.  TP))  THEN 
NMP^NMP*! 

KX(NMP)^KX.X 
KY(KMP;-XYT 
ELSE 
END  !.'■ 

CCNTINUr 

CONTINUE 


2»» 

:9« 

je;  C 

3«H 

390 

300  C 

301  C 

302  C 

303 

304  C 
005  C 
3oe  C 
so:' 

300 

3<I9 

310 

,11 

312  1 

313 

314 

315  C 

316 

317 

310  C 
110  C 
.110  C 

311 
372 
333 
324 
32.1 
320 
327  C 
330  C 
320  C 
,1.10 
3.11  C 

332  C 

333  C 
3.14 
135  C 
.130  C 
.1.17  C 
.135  C 
.lie 

3  I'l 
3  ll 

.112  100 
343  C 


RETURN 

END 

SUBROUTINE  TAPFR(NE.U,T ,X,T ,ET,N .HL.S) 

REAL  0(NE),V(HE),X(»E),T(tE) 

•■••■••or  *•«**«*•«••*  3*  «•*•*«**•* 

*  COMPUTE  APPROXIMATE  ARAAT  RADIUS .  * 


ArS*Nt 


♦  COMPUTE  AMPLITUDE .  * 


DO  1  K=1.»E 

R-SQKT(1(K)*X(K)+T{X)*I(K)) 
AMP-rx+Cl  ,0-ET)*((l  .0-(R/  A>  +  *2)i*»N) 
U(K)=AMP 
V(K)-^0.0 
CONTINUE 
flETORN 
END 


SlTBRCOTINf.  CGRAt)(BE,BSP,NM?,NS,l!l,C,V.X,T,SU.SV,MP,IP,KX,K7,GU,GP 
A  ,aUU,aUV,HVU,UVV,PU,PV,QU,QV,BU,RV,ZU,ZV; 

COMPUTE  POWER  IN  A  SINGLE  DIPECTIOK.  • 


REAL*4  U(N.T)  ,V(tlE)  ,I(NE)  ,Y(NE)  ,SO(NB)  ,SV(ME>  ,OD(NE)  ,CV(NE) 
REAL*'4  HW/(NE,NE),aUV(NE,NE),HVO(NE,»E),HVV(»E,»E) 

REAL*4  PU(KE),rV(ME),yU(NE),QV(llE),RO(NE),RU(BE),ZU(SE),ZV(KE) 
REAL*4  KI(NSP)  ,KT(NSP) 

REAL*4  KXI.KTI 
INTEGER  IP(NSP) 


•  INITIALIZE  THE  NUMBER  OF  SIDELOBE  SAMPLE  POINTS. 


NP  =  0 


•  PERFORM  :.'S  RE-LTARTS. 


* 


DO  1  L=1,NS 


*  FIND  THE  PEAK  SIDELOBE  POWER  POINT  Ai. .  UPDATE  THE  HUMBER  • 

•  OF  MATCH  POINTS. 


CALL  rEAKSL(NE,U,V,I,I,NSP,NMP,KI.KT,IP,NP,PMAX) 
WRITK(*,160}  (IP(I),I=1,HP) 

WRITE(2,lfiO) 

FORKATdOIB) 


•18 


344  C 
34&  C 

346 

347 

345 

349  C 
360  C 
351  C 
353  C 

353  C 

354 

355 

350 

357 

358 
350 

300 

301  2 
363  C 

303  C 

304  C 

305 
300 

307 

308 

369  3 

370 

371 
37J 

373  200 

374  t 

34  5  C 

376  C 

377  C 
370 
3*/0  C 

300  w 

301  C 

302 
3  43 
301 
385 
300 
307 
300 

300  S 
3UI> 

301 

302  C 


FIND  THE  POWER  IK  THB  NilK  BEAM. 


KXI=0.0 

KTI=0.0 

CALL  POWI(NE,B.V,I»T,KII.KTI,PO) 


*  COMPUTE  AVD  STORE  AK  ARIUT  OF  THE  1ST  DERIVATIVE  OF  THB 

«  POWER  IK  THE  MAIK  BEAN  DIRECTIOK  WITH  RESPECT  TO  THE  MTE 

*  VARIABLE. 


KXI=0.0 
KTI=0.0 
DO  2  M=1.KE 

CALL  DPOWIU(NE,0,T,I,T,M,KXI,KT1,PPWIO) 
SU(M)=DPW10 

CALL  DPOWIV(NE.U.V,X.T.M.KXI.KTI.DPWIV) 
SV(M)=DPWIV 
COKTINUE 


*  COMPUTE  AND  STORE  THE  AVERAGE  SIDELOBE  POUER.  * 


PA  =  0.0 

DO  3  1=1, NP 

CALL  P0WI(»E,U,V,X,T,KX(IP(I)),KT(IP(1)),P0WER1) 
PA=PA»P0UER1 
CONTINUE 
PA=PA/NP 

WRITEC* , 200)  L , 10 . 0*ALOG10(PA/P0) , 10 . 0*AL0010( PMAX/PO) 
WRITE(2 , 200)  L , 10 . 0*ALOG10(PA/P0) , 10 . 0»AL0010(PHA1/P0) 
FORMATdI.IB.lX.'AVE  PEAK  SLL' ,lX,FlB.«,iX, 

'PEAK  SLL’ ,lX,F16.e) 


•  COMPUTE  THE  NEGATIVE  OF  THE  GRADIENT. 


DO  4  M=1 ,NE 


•  COMPUTE  THE  DERIVATIVE  OF  THE  AVERAGE  SIDELOBE  POWER.  • 


DPAMU=0,0 
DPAMV=0.0 
DO  6  1=1, NP 

CALL  DFOWrU(«E,U,V,X,T,M,KX(IP(I)) ,KT(IP(I)),DPWIO) 
DPAMU=CPAHU*DPW1U 

CALL  DP0WIV(NE,U,V,X,T,M.KX{IP(I)).KI{1P(I)) ,DPUIV) 
Df AMV=DPAMV4DPWIV 
CO.NTrNUE 
DPA«tI=DPAMy/!fr 
DPAMVrl)FAMV/NP 
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393  C 

394  C 

396 
39« 

397  4 

398  C 

399  C 

400  C 

401 
40] 

403 

404 
406 

406 

407 

408 

409 

410 

41 1 

412 

413 

414 
416 
416  8 
41  7 

418 

419 

420 

421 

422 

423  7 

424  6 

426  C 
428  C 

427  C 

428 

429 

430 

431 

432 

433 

434 
43.6 
438 
437 

4.16  1  1 
419 
440 
-.4, 


*  CONPDTE  TEE  lEGiTITB  GRiDIEHT. 


4 


CU(H)  =  (PA*SU(M)/PO-t)PAME)/PO 
CV(M)=(PA*Sy(M)/PO-DPAMV)/PO 
COSTIHOE 


*  COMPETE  THE  STMKBTRIC  BESSIAI  SDB  MATRICES  ET7U  AID  HTT  * 


DO  0  M=1,NE 
DO  7  I=H,RE 

KXI=0.0  _  _  _  _ 

KT1=0.0 

CALL  ODPltnj(RP.,X,T.M.S.KIl,KTI,DDPDD) 

DDPOaU=DOPOT 

CALL  DDP1VV(HE,X.T,M,S.KII,KTI,DDPVV) 

DDPOVV=DDPVV 
DDPAUU-O.O 
DDPA7V=0.0 
DO  8  1=1, HP 

CALL  DDPHni(HE,X,T.M.S,KI(IP(I)),KT(lP(I)),DDPOO) 
DDPAinj=DDPADD+DDPUU 

CALL  DDPIVV(HE.X.T.M.lt.KX(IP(I)),KT(IP(I)),DDPVV) 
DDPAVy=DDPAVV4DDPVy 
CONTIROE 

DDPAmj=DDPAmJ/RP 

DDPA7V=DDPAVV/RP 

Htnj(H,N)=(CD(K)*SU(M)+GO(M)*SD(R)+DDPAD0-PA*DDPOinj/PO)/PO 

HUD(R,M)=HOU(M,N) 

H7V(M,N)=(GV(R)*S7(M)+G7(M)*S7(H)4DDPA77-PA*DDP077/P0)/P0 

HVV(N,M)=H77(M,H) 

CORTIREE 

CONTINEF. 


*  COMPETE  THE  ASTMMETRIC  HESSIAR  SEB  MATRICE  HE7 .  * 


DO  e  M=1,RE 
DO  10  R=1,HE 
KII=0,0 
K!I=0.0 

CALL  DDPIEV(RE,X,T,M,R,KXI,KTI,DDPE7) 

DDPODV=DDPUV 
DDPADV=0.0 
DO  11  1=1, RP 

CALL  DDPIOV(RE,I,T,H,B,XX(IP(I)).KT(IP(I)),DDPOV) 
DDPAU7=DDPAUV4DDPUV 
CONTINUE 

DDPAUV=DDPAUV/NP 

HUV(M,N)=(OV(N)*SU(M)+GU{M)*S7(N)+DDPAE7-PA*DDP007/P0)/P0 

HVU(N,M)=HUV(M,N) 
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442  10 

443  0 

444  C 
«45  C 
44«  C 

447  C 

448  C 

449  C 

450 
481 
452 

483 
454 

465  12 

458  C 

467  C 

468  C 

459 

480 

481 

482 
463 

484 

485 

466 

487  13 

488  C 

469  C 

470  C 

471 

472  C 

473  C 

474  C 

476 
4  76 

477 

478 

479 
180 

481 

482 

483 

484 

485  IB 

486 

487 

488 

489 
400 


CONTIHUE 

CORTIHUE 


*  START  THE  CONJUGATE  GRADIENT  ALGORITHM. 

* 

*  INITIALIZE  TEE  RESIDUAL. 

* 

DO  12  H=1,NE 

RU(M)=GD(M) 

RV{M)=GV(M) 

ZU(M)=0.0 

ZV(M)=0.0 

CONTINUE 

*  INITIALIZE  SEARCH  VECTOR. 

* 

CALL  MATVEC(RE.HUU,R0,HDV,R7,Q0) 

CALL  MATVEC(NE,HVO,RU,HVV,RV,QV) 

CALL  N0RM22(RE,Q0,qUN) 

CALL  N0RM22(NE,qV,qVN) 

BE0=1 .0/(QUN+QVN) 

DO  13  1=1, HE 

PU(I)=BEO*qO(I) 

PV(I)=BEO*qV(I) 

CONTINUE 

*  PERFORM  CONJUGATE  GRADIENT  ITERATIONS. 

* 

DO  14  K=1,NI 

*  UPDATE  AMPLITUDE  VECTOR  AND  RESIDUAL. 

* 

CALL  MATVEC(HE,mnj,PD,HOV,PV,qU) 
CALL  MATVEC(NE,HVO,PU,HVV.PV,QV) 
CALL  N0RM22(RE.QU,QUN) 

CALL  R0HM22(RE,qV,qVR) 

AK=1  .0/(qUN4-qVN) 

DO  16  1  =  1, RE 

ZU(I)=ZU(I)+AK*PD(I) 

ZV(I)=ZV(I)+AK*PV(I) 

RU(I)=RU(I)-AK»qrj(I) 

RV(I)=RV(I)-AK*qV(I) 

CONTINUE 

CALL  N0RM22(NE,RU,RDNM) 

CALL  N0RM22(NE,RV,R7NM) 
ERR-SQRT(RDRM+RVNM) 

WRITF,(*,600)  L,K,ERR 
WR1TE(2,B00)  L,K,ERR 
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4»1  COO  FOC.UTOX. 'KBSIDOiL' .1X,I3.1X,I3.1X.B16.8) 

492  C  ••***•*****»•,,*,,**»•,*••***••♦♦•**••♦**•***♦•***•*****•*•*** 

403  C  *  UPDiTE  SEARCH  VECTOR.  * 

494  C  ♦♦*4*»**4i********#*9*********9***»****99***9***********9****** 

495  CALL  NATVEC(RE,E01I,RO,HUT,RV.qD) 

496  CALL  NATVEC(NE.HV0,R0,HV7.RT,qV} 

497  CALL  H0RM22(VE,q<J,QUH) 

498  CALL  l(0RM22(EE.qV.qVN) 

489  BEK=1.0/(qUN+qVB) 

500  DO  16  1=1, HE 

501  PU(I)=P0(I)4-BEK*QD(I) 

502  PV(I)=PV(I)+BEK*QV(I) 

503  16  CONTIHOE 

504  14  CONTINUE 

505  C  •**••»•*♦*••*♦*♦*********♦*•*•***•************♦•****♦**•••♦*»* 

508  C  *  UPDATE  AMPLITUDE  VECTOR.  • 

507  C  ♦**♦♦♦»**♦*****♦*♦♦♦***••♦*♦*♦♦******♦*********•**•♦*♦***♦**•* 

508  DO  17  1=1. HE 

509  U(I)=D(I)+ZO(I) 

610  V(I)=V(I)4ZV(I) 

611  17  CONTINUE 

612  CALL  NORMAL(NE,U,V) 

313  C  ***•♦**•**********♦***•*•********•**•******♦**•**♦♦**♦*♦••**** 

514  C  *  WRITE  OUT  THE  NEW  ESTIMATE  TO  A  FILE.  • 

515  C  **•*••****••••••*•*•*•••••*••*•*•*******•»••♦********•****••*• 

510  REWIND  4 

617  WRITE(4,101)  (I,0(I),V(1).I=1,NE) 

518  101  F0RMAT(IB,1X,F1B.7,1X,F16.7) 

619  IF  (NP  ,EQ.  NMP)  THEN 

520  GO  TO  00 

621  ELSE 

622  END  IF 

523  1  CONTINUE 

524  09  CONTINUE 

625  RETURN 

528  END 

627  C 

628  SUBROUTINE  PEAKSL(NE .0, V ,X,T .NSP.NMP.KX.KT, IP.NP ,PMAX) 

629  C  ♦***•♦♦**♦•****•*****•*♦•**♦•******♦♦**»**♦**♦**♦•*•***♦**♦•••*••** 

530  C  *  THIS  SUBROUTINE  FINDS  THE  PEAK  SIDELOBE  POWER  LEVEL  COORDINATE  * 

531  C  *  POINT  (KX(LMAX),KT(LMAX)).  IF  THE  POINT  IS  HEW  THEN  IT  IS  * 

532  C  •  ACCUMULATED  INTO  THE  POINTER  ARRAY  IP.  THE  HUMBER  OF  POINTS  IN  * 

633  C  *  THE  INTEGER  POINTER  ARRAY  IP  IS  NP  SUCH  THAT  IP(I) ,1=1. . .NP.  * 

634  C  *»****•***»♦•••♦*•****•***♦♦*♦**•***♦**•*••»*♦*♦***♦*♦*»*♦***♦*•**♦ 

635  REAL  U(BZ),V(NE),X(NE),T(NE),KX(NSP),KY(HSP) 

638  INTEGER  IP(NSP) 

537  C  *••♦♦*♦••*♦♦»**•**••*•••♦***•***•*****•***«••**•»***«****•**** 

638  C  *  COMPUTE  THE  POWER  AT  EACH  POINT  AND  FIND  THE  MAX.  * 

539  C  ♦♦****»********•*•***♦♦***♦•**♦♦*•»**********•*♦********♦***•* 
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640 


641 

642 

643 
544 
546 
646 
847 
846 
640  1 
660  C 
BSl  C 

552  C 

553  C 

554  C 

655  C 

656 

657 
556 

550 
560 

set 

6B2  2 

563 

664 

sds 

666 

507 

566 

660 

670  C 

571 

672  C 
573  C 
674  C 

576 
570 

577 
670 
570 
SHl) 

SBI 

5SJ 

5R3 

551  1 

555 

556 

557 

558 


PMAX=-3000.0 
LMiX=-l 
DO  1  1=1, RMP 

CALL  POWI(HE,D,V,I,T,KX(I) ,KT(I) ,PW1) 
IF  (PWI  .GT.  PMAI)  THEM 
LMAX=I 
PMAX=PWI 
ELSE 
END  IF 
CONTIMHE 


*  COMPARE  LMAI  TO  ALL  THE  PREVIODS  POISTERS  IP(1) , 1=1 . . .NP . 

•  IF  LMAX  MATCHES  OHE  OF  THE  PREVIOOS  POIHTERS  THEN  THE 

•  POINTER  ARRAY  HEED  NOT  BE  DPDATED  (1FLAG=1) .  IF  LMAX  DOES 

*  NOT  HATCH  A  PREVIOUS  POINTER  VALUE  (IFLAG=0)  THEN  IT 


IFLAC=0 
DO  2  1=1, NP 

IF  (IP(I)  .EQ.  LMAX)  THEN 
IFLAC=1 
ELSE 
END  IF 
CONTINUE 

IF  (IFLAG  .Eg.  0)  THEN 
HP=HP+1 
IP(NP)=LMAI 
ELSE 
END  IF 
RETURN 
END 

SUBROUTINE  POWI(NE .0 ,V ,I ,T,KII ,KTI .POWERI) 


*  COMPUTE  POWER  IN  A  SINGLE  DIRECTION.  * 


REAL  U(NE) ,V(NE).X(NE).T(NE).KXI,KTI 
SUM1=0.0 
SDM2=0.0 
DO  1  K=1 ,NE 
PSI=KXI»I(K)+KTI*T(K) 

CP=C0S(PS1) 

SP=SIN(PSI) 

S0M1=SUM1+D(K)*CP-V(K)*SP 

SDM2=SUM2+U(K)*SP+V(K)»CP 

CONTINUE 

TIrSDMl*SUMl 

VI=SUH2*SUM2 

POWERI=TI+VI 

RETURN 
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SS9 

seo  C 

SPl 

(93  C 
893  C 
(98  C 
(98  C 
896 
807  C 
BOB  C 
899  C 
BOO 
BO  I 
BO  3 
eo3 
BO  4 
605 
BOB 
BO  7 
BOB  1 
BOB  C 
BlO  C 
flu  C 

B13 
Bi3 
BM 
615 
BIB  C 
B17 
BIB  C 
619  C 
(130  C 
621  C 
022 
623  C 
634  C 

626  C 
626 

627 
fl2» 

629 

630 

631 

632 

633 

634  1 

635  C 

636  C 

637  C 


END 

SUBROUTINE  DPOUIU(NE ,0 , ? , 1 ,T . M .KII ,KTI .DPWIO) 

•  COMPUTE  1ST  DERIVATIVE  OF  POWER  WITH  RESPECT  TO  THE  MTH  VARIABLE* 

•  IN  THE  ITH  DIRECTION.  * 

REAL  U(KE) , V(HE) ,X(HE) ,T(HE) .KII.KTI 

•  COMPUTE  SUMS.  * 


SDM1=0 .0 
SUM2=0 .0 
DO  1  K=1,NE 
PS1=KII*X(K)+KTI*T(K) 
CP=COS(PSI) 

SP=^S1B(PSI) 

SDM1=SUM1>D(K)*CP-V(K)*SP 

SUM2=SUM2*U(K)«SP+V(K)*CP 

CONTINUE 


•  COMPUTE  1ST  DERIVATIVE  OF  POWER.  * 


PSIM=KII*I(M)4-KTI*T(M) 

DPV,’IU=2.0*(COS(PSIM)*SDM14SIN(PSIM)*SOM2) 

RETURN 

END 

SUBROUTINE  DPOWIV(NE,U,V,I,T,M,KXI,KTI,DPWIV) 


•  COMPOTE  1ST  DERIVATIVE  OF  POWER  WITH  RESPECT  TO  THE  MTH  VARIABLE* 

•  IN  THE  ITH  DIRECTION.  * 


REAL  U(NE) ,V(NE) ,X(HE) ,T(NE) ,KXI ,KTI 


*  COMPUTE  SUMS. 


« 


S0M1=0.0 
SDM2=0,0 
DO  1  K=1,NE 
PSI=KX1*X(K)+KTI*T(K) 
CP=COS(PSI) 

SP=SIIUPSI) 

SUM1=SUM1+U(K)*CP-V(K)*SP 

SUM2=SUM2+U(K)*SP+V(K)*CP 

CONTINUE 


•  COMPUTE  1ST  DERIVATIVE  OF  POWER. 


* 
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03S 
«39 
«40 
341 
84J  C 

343 

344  C 

345  C 
343  C 
347  C 
348_  . 
349  C 
330  C 

351  C 

352 
853 
054 
055  C 
656  C 
857  C 
058 
859 

830 

831  C 
062 

303  C 

304  C 
035  C 

630  C 

007 
388  C 
009  C 

670  C 

671 
072 
673 
074  C 
075  C 
070  C 
077 
878 
070 

080  C 
681 

002  C 

003  C 
094  C 

006  C 
000 


PSIM=KXI*I(M)4KTI*T(M) 

DPUIV=2.0*(C0S(PSlM)*SnM2-SIH(PSIM)*SDm) 

RETURN 

END 

SUBROUTINE  DDPIOU(HE ,X .T , H . B ,KII ,KTI .DDPUU) 


•  COMPUTE  2ND  DERTVXTITE  OF  POWER  WITH  RESPECT  TO  THE  M  AND  NTH  * 

•  VARIABLES  IN  THE  I^:H  DIRECTION.  • 


REAL4-4  X(NE),T(NB),KX1.KTI 


0  COMPUTE  DIFFERENCE  OF  PSI  FUNCTIONS. 


* 


PSIM=KXI*I(M)4-KTI*T(M) 

PSIN=XXI*X(NUKTI*T(N) 

DIF=PSIM-PSIN 


*  COMPUTE  2ND  DERIVATIVE  OF  POWER. 


* 


DDPDU=2.O0COS(DIF) 

RETURN 

END 

SUBROUTINE  DDPIUV( NE , X .T , M ,N ,KXI ,KTI .DDPUV ) 


•  COMPOTE  2ND  DERIVATIVE  OF  POWER  WITH  RESPECT  TO  THE  M  AND  NTH  * 

•  VARIABLES  IN  THE  ITH  DIRECTION.  • 


REAL*4  X(HE),T(HE),KXI,KTI 


•  COMPOTE  DIFFERENCE  OF  PSI  FUNCTIONS. 


* 


PSIM=KXI*I(M)+KTI*T(M) 

PSIN=KXI*X(N)+KTI*T(N) 

DIF=PSIM-PSIN 


*  COMPUTE  2ND  DERIVATIVE  OF  POWER. 


* 


DDPDV=2.0»SIN(DIF) 

RETURN 

END 

SUBROUTINE  DDPI VV ( NE , X , T . M .N .KXI ,KTI , DDPVV ) 


0  COMPOTE  2ND  DERIVATIVE  OF  POWER  WITH  RESPECT  TO  THE  M  AND  NTH  * 
•  VARIABLES  IN  THE  ITH  DIRECTION.  • 


REAL04  X(HE)  ,T(NE),KXI,KTI 
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#87  C  *•**»••**•*••*••**♦**•****•••**•*•••*••♦•*•*•*♦****•*•*♦**••*••••*• 

ess  C  *  COMPOTE  DIFPBREHCE  OF  PSl  FOSCTIOHS.  • 

ese  C  ♦*•*•?*•**♦**•****♦•*****•*****♦•♦♦ 

ewo  PStM=KXI*X(M)4KTI*T(M) 

e#i  PSIH=KXI*X(K)+KTI*T(H) 

092  DIF=PS1M-PSIB 

093  C  *••**•♦♦*•»♦•****••***••*********♦ 

094  C  *  COMPOTE  2BD  DERlViTTVE  OF  POWER. 

090  C 

090  DDPTV=2.0*C0S(D1P) 

097  RETORN  .  _  _ 

098  END 

099  C 

700  SOBROOTINE  H0RM22(R,i.iN) 

701  C  **•**••**•••••**•*•#•*••*•**•***•• 

702  C  *  THIS  SOBROOTINE  COMPUTES  THE  EOCLIDIIH  NORM  SQOARED  OF  A.  ♦ 

703  C  ♦»»**•*••«*••••******••**•••••••**•*••*•••*••****•••***•****•••*••* 

704  REAL*4  A(N) 

705  AN=0.0 

706  DO  1  1  =  1,  N 

707  Ai;  =  AN+A(I)*A(I) 

700  1  CONTINUE 

709  RETURN 

710  END 

711  C 

712  SOBROOTINE  MATVEC(N,A,B,C,D,E) 

713  C 

714  C  •  THIS  SUBROUTINE  COMPOTES  THE  MATRIX  PRODOCT  E=AB+CD.  • 

715  C  •**•***••*•*******•****•••••••*••***•*••**•••**•*•**••••••***••*••• 

716  REAL*4  A(H.N).B(H).C(N,N).D(N).E(N) 

717  DO  1  1  =  1, N 

718  SUK=0.0 

719  DO  2  J  =  1,N 

720  SOM=SUM+A(I,J)*B(J)+C(I,J)»D(J) 

721  2  CONTINUE 

722  E(I)=SOM 

723  1  CONTINUE 

724  RETURN 

720  END 

726  C 

727  SUBROUTINE  NORMAL(N ,A,B) 

728  C  *••**•••*•***•••**••«***•*••**••*••*••*••••***••*»*****♦*•••♦•»•**» 

729  C  *  THIS  SOBROOTINE  NORMALIZES  A  VECTOR  OF  LENGTH  N.  * 

730  C 

731  REAL*4  A(N),B(N) 

732  VMAX=50RT(A(1)*A(1 , *B(1)*B(1)) 

733  DO  1  1  =  1, N 

734  V=SQRT(A(I)*A(I)+B(1)*B(I)) 

735  IF  (V  .GT.  VMAX)  THEN 
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736 

737 


738 

739  1 

740 

741 
74J 

743  2 

744 

745 

746  C 

747 
74B  C 
740  C 
760  C 
751 
,f2 

753 

754 

755  C 

756  C 

757  C 
768 
750 

760 

761 

762 

763 

764 

765 

766  C 

767  C 

768  C 
760 

770  C 

771  C 

772  C 

773 

774 

775 

776 

777 

778 
770 

780 

781  2 

782  1 

783 

784 


VMAI=V 
ELSE 
END  IF 
CONTINUE 
DO  2  1=1, N 
A(1)=A(I)/ABS(VMAI) 
B(I)=B(I)/ABS{VMAX) 

CONTINUE 

RETURN 

END 

SUBROUTINE  CAIN(NE.0,7.1,T,GDB) 


*  COMPUTE  POWER  IN  A  SINGLE  DIRECTION.  * 


REAL  U(HE),V(NE).X(NE).I(KE) 
RAD=.17463293E-01 
PI= .3141B93E+01 
TP  .02e318SE+01 


•  GRAPH  INPUTS  * 


NP=41 

NR=31 

PMIN=0.0 

PMAX=360.0 

TMIN=O.C 

TMAX=90.0 

DP=(PMAX-PM!N)/(NP-1) 

DR=(T«AX-TMIN)/(RR-1) 


*  COMPUTE  PEAK  VALUE. 


* 


CALL  POWER(NE ,U ,V ,X ,T,0 . 0 ,0 . 0 .PEAK) 


•  COMPUTE  THE  SOLID  ANGLE. 


SUM=0.0 

DO  1  I=0,NR-2 

T--T?1IN+DR/2.0  +  I4'DR 
ST=SIN(RAD*T) 

DO  2  J=0,NP-2 

P^PMIN^DP/2.0+J*DP 
CALL  POWER(NE,U,V,X,T,T,P.POW) 
snM=snM+ST»pow 
CONTINUE 
CONTINUE 

TEST=RAD*RAD*DR*DP*SUM 

S0L1D=RAD*RAD*DR*DP*S0M/PEAK 
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Tse  GDB=10.0*iL0010(TP/SOLlD) 

786  RETORH 

787  END 

788  C 

78B  SOBROUTINE  POWPAT(IE.O,V ,1,T) 

7»0  C 

7B1  C  *  COMPUTE  POWER  IR  k  SIRGLS  DIRECTION.  • 

7»j  C  •••**•***••*•*••••*•*•••••*•• 

7B3  REAL  D(HE).V(HE).I(RE).T(RE) 

791  OPEB(OKIT=7.FILE=’DATl*) 

796  RAD=.174B32B3E-01 

798  P1=.3141B93E+01 

797  TP= .fl283186E+01 

798  C  **♦******•**•****•*♦*♦♦*•♦*•* 

799  C  ♦  GRAPH  INPUTS 

800  C  •****•••*»••••••*••*••*••••*• 

801  NP=6B 

802  NR=30 

803  KS=0 

804  pmN=o.o 

80P  PMAI=360.0 

806  TMIN  =  0.0 

807  TMAX=90.0 

808  DP=(PMAI-PMIN)/(NP-1) 

808  DR=(TMAX  ■TMI5)/(NR-1) 

8;c  AI-160.0 

811  AT=AI 

812  AZ=fl80.0 

813  PO  =  BB.O 

814  Q0=60.0 

815  IA  =  -2.0 

816  IB-2.0 

817  TA  =  -2.0 

818  7B=2.0 

810  ZA  =  -2 . 0 

820  ZB -2. 3 

821  AT=.3 

822  WRITE(7,*)  NP.NR.NS.AI.AT.AZ 

623  WRITE(7,*)  XA,XB,TA,TB,ZA,ZB 

824  WRITE(7,*)  AT,PO,qO 

825  C  ♦•♦••*»****««***«************i 

828  C  *  PHI  VALUES.  * 

827  C  •.••*»»«*«»«*»»**»«*»»*»»***»»»*******»*»*»*****»**»*i*»************ 

828  DO  1  I=0,NP-1 

820  P=PMIN-*I*DP 

830  WRITE(7,»)  RAD'P 

831  1  CONTINUE 

832  C  •*•♦•*••*•»*•♦*••*»*•*♦*•••••*•*••*•*****••**••**•••*•••••*•••**••* 

833  C  *  RADIAL  (THETA)  VALUES.  • 
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S3S 

DO  2  1=0,IJR-1 

T=TMIH+1*DR 

ear 

WRITE(7,*)  RAD*T 

eas  2 

CONTIBUE 

830  C 

840  C 

*  FDBCTION  VALUES. 

841  C 

843  C 

843  C 

*  COMPUTE  PEAK  VALUE. 

844  C  . 

946 

CALL  P0WER(HE,U,V,I,T,O.O,O.O,PEAX) 

848  C 

847  C 

*  COMPUTE  BACH  VALUE. 

* 

848  C 

040 

DO  3  I=0,NR-1 

0SO 

T=TMIN+I*DR 

851 

DO  4  J=0,BP-1 

053 

P=PMIB+J*DP 

053 

CALL  P0WER(HE,O,V,I.T,T.P.P0W) 

061 

WRITE(7,*)  POW/PEAK 

855  4 

CONTINOE 

858  3 

COSTIBUE 

857 

RETURB 

858 

EBD 

859  C 

880 

SUBROUTINE  SPOWPAT(NE,U, V.I .T) 

881  C 

883  C 

•  COMPOTE  POWER  IN  A  SINGLE  DIRECTION. 

* 

883  C 

004 

REAL  D(NE) ,T(BE),X(BE) ,T(HE) 

SflS 

0FEN(DNIT=8,FILE-'SDATA’) 

600 

RAD=.174B3293E-01 

007 

PI=.3141683E+01 

866 

TP=.8283186E+01 

889  C 

870  C 

•  GRAPH  INPUTS 

0 

871  C 

^^0009 

872 

BP=11B 

873 

NR=63 

874 

NS=0 

875 

PMIN-0,0 

870 

PHAX--3flO  ,0 

877 

THIN=0.0 

878 

TMAX^eo.O 

879 

DP=(PMAX  PMIN)/(NP-1) 

8#0 

DR=(TMAX-THIN)/(NR-1) 

80! 

AX=800.0 

802 

AT-AX 
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883 

iZ=AZ 

684 

PO=2B.O 

88S 

00=60.0 

886 

Ii=-2.0 

88’ 

IB=2.0 

888 

o 

(N 

II 

88V 

TB=2.0 

890 

ZA=-2.0 

891 

ZB=2.3 

893 

AT=.3 

883 

IT=4 

894 

WRITE(8.*)  NP,HR.HS,IT.4I.1T,AZ 

605 

WRITECe.*)  IA,IB.TA,TB,ZA.ZB 

886 

WRITE(8.*)  AT.PO.qO 

887  C 

888  C 

*  PHI  VALUES. 

9 

6B8  C 

900 

DO  1  I  =  0,NP--1 

901 

P=PMIN+I*DP 

003 

WRITE(e,*)  RAD*P 

9C3  1 

CONTINUE 

e04  C 

BOS  C 

‘  RADIAL  (THETA)  VALDES. 

* 

ene  C 

907 

DO  2  I=0,NR-1 

908 

T=TMIN+I‘DR 

900 

WRITE(8,*)  RAD»T 

BIO  2 

COMTINOE 

ei  1  C 

B13  C 

•  FUNCTION  VALDES. 

* 

813  C 

911  C 

BIS  C 

♦  COMPUTE  PEAK  VALUE. 

* 

916  C 

oir 

CALL  POWER(NE.U,V,I.T. 0.0, 0.0, PEAK) 

918  C 

BIB  C 

•  COMPUTE  EACH  VALUE. 

* 

020  C 

921 

DO  3  I=0,NR-1 

9?r 

T=TMIN+I»DR 

9  23 

DO  4  J=0,NP-1 

924 

P^PMIN-fJ^DP 

925 

CALL  POUER(NE,U,V,I,T,T,P,POW) 

020 

VRlTECe,*)  POW/PEAK 

B27  4 

CONTINUE 

B28  3 

CONTINUE 

u:i’. 

RETURN 

93f 

93!  C 

END 

6fi 


032 
033  C 
03-*  C 
035  C 
030 
037 
038 
030 
040 

or. 

043 

043 

'<>44 

0 

048 
Oil’ 
051,- 
OBI  1 
or: 

053 

S^4 

ers 

050 


SOBRO0TINE  POWER( HE ,0,V,X,T,T,P,POW) 


«  COMPtTTE  VOUER  IH  i  StHGLE  DIRECTIOH  GIPEH  i  THETA  AND  PHI.  * 


REAL  U(NE) ,V(HE),X(HE).T(HE),K1I,KTI 

RAC^.174E32C3B-0i 

TP=.e2e3106E+0l 

TPS=TP*SIH(RAD*T) 

RP-RAD*F 

KXI=TPS*COS(RP) 

KTI=TPS*SIH(KP) 

SDMI -0.0 
SOM2-0.0 
DO  1  K=1,NE 
PSI=KXI*X(K)tKTT*T(K) 

CP-COS(PSjC) 

sp  ^sifjrpsi) 

SUMl^SUfll+U(K)*CP-V{K)*SP 
aoM2  =  SUH2+0(K)*SP  +  V(K)'‘CP 
COHTIHUE 
TI=SUM1»3UM1 
VI=S0M2*SDM2 
PO«’-TI<V,: 

RETURi; 

END 
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PROGRAM  DIFSTH 


5  C 

3  C 

4  C 
t  C 
0  C 

7  C 

3  C 

p  c 

10  C 

11  c 
17  C 
13  C 
I  4  C 

15  C 

16  C 

17  C 
IS  C 
10  C 

20  C 

21  C 

22  C 

23  C 

24  C 

25  C 
20  C 

27  C 

28  C 
20  C 

30  C 

31  C 

32  C 

33  C 
31  C 

35  C 

36 

37 

38 
30 
40 
4  1 

42 

43  C 

44  C 

45  C 
4fl 


47 


40 


THIS  PROGRAM  ST5TBESIZES  THE  AMPLITUDE  AHD  PHASE  DISTRIBDTIOH 
NECESSiRT  TO  GENERATE  A  DIFFERENCE  PATTERN  07ER  A  HEXAGONAL 
mAT  OF  POINT  SOURCES  WITH  SOME  ELEMENT  FAILURES. 


TlMOTilT  J.  PETERS  ' ‘ST  UPDATED 

THE  AEROSPACE  CORPORATION  ./Bl 

EAST  EL  SBGUNDO  BOULEVARD. 

EL  SEGUNDO,  CA  B024B 


INPUTS: 

NL  -  NUMBER  OF  LAYERS  FORMING  THE  BEX AGONAL  LATTICE 

NTE  -  HUMBER  OF  TOTAL  ELEMENTS  =(l+3*NL*(NL+l)-NL)/2  FOR  NL  ODD 

=(3*»L*(HL41)-KL)/2  FOR  NL  EVEN 

NBE  -  NUMBER  OF  BAD  ELEMENTS 

BAD(1..NBE)  -  INTEGER  ARRAY  HOLDING  THE  NUMBER  POSITION  OF  EACH 
BAD  ELEMENT  IN  TEE  ARRAY. 

DS  -  THE  ELEMENT  SPACING  IN  HATELENGTBS. 

NE  -  NUMBER  OF  ACTUAL  ELEMENTS  =  H7E-NBE. 

SLL  SIDELOBE  LEVEL  OP  DESIGN  BATLISS  DISTRIBUTION. 

NBAR-  aRAMETER  FOR  BATLISS  DISTRIBUTION. 

BBF  -  THE  BEAM  WIDTH  BROADENING  FACTOR.  THIS  ALLOWS  THE  BEAM 

TO  BROADEN  BEYOND  THE  DESIGN  VALUE  IN  ORDER  TO  REDUCE  THE 
SIDELOBES  WHE;;  ELEMENTS  FAIL. 

NSP  -  NUMBER  OF  SAMPLE  POINTS  IS  THE  SIDELOBE  REGION. 

MS  -  NUMBER  OF  RESTARTS. 

NI  -  NUMBER  OF  CONJUGATE  GRADIENT  ITERATIONS  PER  RESTART. 
OUTPUTS : 

U(NE)  -  REAL  PART  OF  THE  EXCITATION  OF  EACH  ELEMENT. 

V(NE)  -  IMAGINARY  PART  OF  THE  EXCITATION  OF  EACH  ELEMENT. 


PARAMETER  (ML=6 ,STE=43 , NBE=0 ,NE=43 , NSP=700) 

REAL»'l  U(NE).X(NE),T(NE).S(NE).CH(KE),H(KE,NF) ,P(NE),QC»E) .R(KE) 
REAL*4  Z(KE).KX(NSr),KT(NSP),KIP,XTP,NWBP 
INTEGER  IP(NSP)  ,BAD(NE) 

PI=  ,3141B93E401 

CPEn(UNIT=2,FILE= 'OUTPUT' ,STATUS= 'UNKNOWN' ) 

OPEB(ONIT=4,FILE= 'ESTIMATE' ,STATUS= 'UNKNOWN ' ) 


•  GENERATE  THE  GEOMETRY, 


* 


DS=0.6 

BAD(1)=1S 

BAD(2)^17 

BAP(3)=20 
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so  CiLL  CEONET(HL,DS,RE,HBE,BAD,X,T) 

5,  C  .*•***,♦♦♦.♦*♦♦♦„***♦♦*.♦♦**•*••♦*••*•♦•****.••*♦»♦**.•.•***♦•*♦**• 

5j  C  *  IBITIILIZE  THE  TAPER  DISTRIBOTIOH.  • 

83  C 

ii  IFLACl-l 

66  SLL--40.0 

so  IIBiR=iO 

87  IF  (IFLACl  .EQ.  1)  THEN 

68  CALL  TAPER(KE,D.£,T,SLL,BBAR,IL,»S) 

68  REWIND  4 

80  WRITE(4,101)  (I,D(I),I=1.HE) 

01  -  ELSE  -  -  .  -  _ 

82  END  IF 

83  REWIND  4 

84  READ(4.101)  (I.U(I) ,I=1.HE) 

86  C  **••****♦♦•*•***♦*•**•****•**•**••••••*•*•••••••*••••*•*•♦**••*•••* 

86  C  ♦  COMPUTE  THE  GAIN  (DIRECTIVITT)  OF  THE  IRRAT .  * 

88  C  CALL  GAIN(BE,U,X,T,CDB) 

80  C  **»**•**♦****♦•**•**•♦♦•♦*»*•*•*•••**•••••••*•***•••*•**••*•••••••• 

70  C  *  FIND  THE  HALF  POWER  NULL  WIDTH  BETWEEN  PEAKS.  * 

71  C  . . . . . . . 

72  C  CALL  BEAK(NE,0,X,T,HWBP) 

73  C  . . . 

74  C  *  GENERATE  THE  K-SPACE  SAMPLE  POINTS.  • 

76  C  •••••*•***«»•••••••»*••••• ••••*»**t* 

78  BBF=1.24 

-7  CALL  K3PACE(NMP,NL,DS,SLL,NBAR,BBF.KX,KT) 

78  C  **•**•***••««*••*••**•**•*•«**»***•* ***•**•*****•***«•**•*•******** 

70  C  •  COMPOTE  NORMALIZED  PEAK  SIDELOBE  POWER  IK  DB.  • 

80  C  **••♦*♦•*•♦«*•**•*♦«•«••*»•*••*•**•»•*»****•••*•*•♦***•*•»•»••**•*• 

81  CALL  PEKPOUINE.U.X.T.KXP.KTP.PO) 

82  NP=0 

83  CALL  PEAKSL(NE,D,X,T,NSP,HMP,KX,KT,IP,Hr,PMAX) 

84  rHAXDB.-^10,0*ALOC10(PMAX/PO) 

86  C  •♦••*»♦••*•»••»»••*•»•»••*••*•*••••*•*•»*•******»**•****••••••••**• 

88  C  *  PERFORM  CONJUGATE  GRADIENT  STEPS.  * 

87  C  »•••*•**•** •i********* **♦*•*♦•**••*•* **«8****»********** **♦•••***♦** 

88  NS^IBO 

80  NI=3 

00  CALL  CGRAD(KE,N5P.NMP,NS,NI,U,I,T,S,Nr,IP,KX,KT,GN,H,P,q,R,Z> 

01  REWIND  4 

02  WRITE(4,101)  (I.U(I),I=1,NE) 

0.1  C  4  •••••••*»**•»*»*•***•••»♦•••••***«**»* 

04  C  *  POWER  PATTERN  GRAPHICS  OUTPUT.  * 

95  C  •*»»*»»4*******»«***»*«»***ili*»****4*»t»**t******4*r»»*»**»»4#.*»***« 

08  REWIND  4 

or  READ(4,101)  ( 1 ,0( 1 1  , I=l , NL) 

08  C  •»»»»**»*•** -****«*****0****»***»*»«******4*»*.*»******w«»**«****** 
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•B  C 

*  PRINT  TBB  POWER  PATTERN  IN  CTLIISRZCAL  COORDINATES. 

IGO  C 

lOl 

CALL  P0WPAT(NE,D,X,T) 

to:  C 

103  C 

*  PRINT  TEE  POWER  PATTERN  IN  SPHERICAL  COORDINATES. 

lo-i  C 

:os 

CALL  SPOWPlT(HB,a.X.T) 

)0«  C 

tor  C 

•  FORMATS. 

los  C 

loe  100 

rORMAT(IS,F16.6.1X,Pl6.6) 

-■  - 

no  101 

F0RMAT(I6,P1B.7)  -  -  -  -  - 

........ 

1  1  1 

END 

IIJ  c 

113 

SUBROUTINE  GEOHET ( NL , S , NE .NBE.BAD . X . T ) 

IM  C 

ns  C 

•  COMPUTE  THE  POSITION  OF  EACH  ELEMENT  IN  A  BEXOGONAL  ARRAY. 

na  C 

117 

REAL*4  I(NE),T(NE) 

•  1 

118 

INTEGER  BAD(NE),IFLAC 

■| 

119 

HCT=(SqRT(3.‘))/2.0)*S 

120  C 

121  C 

*  CENTER  hOV. 

* 

122  C 

123 

12  i 

L=0 

1  20 

XMIW=-NL*S 

126 

DO  1  I=0,2*NL 

127 

XX=XMINiI»S 

« 

126 

IF  (XX  .GT.  0.001)  THEN 

12H 

L-L+1 

130 

CALL  CHECK(L,NE,NBE,BAD,IFLAG) 

131 

IF  (IFLAi  .EQ.  0)  THEN 

132 

K=Kn 

133 

I(K)=II 

134 

T(K)r0,0 

loO 

ELSE 

136 

END  IF 

137 

ELSE 

138 

END  IF 

13B  1 

CONTINUE 

Ml  C 

*  TOP  ROWS. 

1  '1  2  V 

1  4  3 

INUM^2»NLtl 

1  -I'l 

DO  2  J=1,UL 

14. *2 

INUM-INUK-l 

• 

1  16 

XMIN^  -m.»SsJ*S/2.0 

1  17 

DO  3  I=0,INUM-l 

■rygVTT’tavj;*?  .  '.vj  r-'v-fry  uj 

G-J 

MS 

xi=iMiir+i*s 

MO 

IF  (XI  .GT.  0.001)  IHEir 

1(0 

L=L+1 

1(1 

CALL  CHECK(L.XE.IBE,BA0,1FLAC) 

163 

IF  (IFLAG  .Sq.  0)  THEH 

163 

K=K+1 

U4 

X(K)=IX 

166 

T(K)r=HCT*J 

15S 

ELSE 

15? 

END  IF 

168 

ELSF, 

169 

END  IF 

leo  3 

CONT.TNUE 

lai  2 

CONTIHOE 

102  C 

too  C 

•  BOTTOM  ROUS. 

1.31  C 

L  03 

IK0>!^2-I}L+1 

leo 

DO  4  J=1,NL 

187 

IBDM-INUM-l 

18^ 

IMIN--NL*S+J*S/2.0 

ISO 

DO  B  I=O.IBOM-1 

1  7n 

ii-iMinn*s 

171 

IF  (IX  .CT.  0.001)  THEN 

172 

1  73 

CALI,  CHECK(L.SE,KBE,EiP,IFLAC) 

I  74 

IF  (IFLAG  .Eg.  0)  TEES 

1  75 

K=K  +  1 

178 

I(K)=II 

177 

T(X)=-HCT»J 

178 

ELSE 

179 

END  IF 

180 

EISE 

181 

END  IF 

182  6 

CONTIROF. 

183  4 

CONTINUE 

184 

RETURN 

186 

END 

186  C 

187 

SUBROUTINE  CHECK( L , HE , HBE ,B AD , IFLAG ) 

188  C 

180  C 

♦  CHECK  TO  SEE  IF  THAT  ELEMENT  IS  BAD. 

9 

190  C 

iw; 

INTEGER  BAD(NE) 

102 

IFLAC^O 

lf(.1 

DO  1  1  1 ,NBF 

(PI 

IF  (L  EQ.  BAD(I))  THEN 

1  or. 

IFI.AG-^l 

i  98 

ELSE 

05 


197  EKO  IF 

198  1  COKTINOE 

199  RETUKN 

300  END 

201  C 

302  SOBRODTINE  BElH(NE,tI,Z,T,BVFN} 

203  C  »••••••*•***•**♦*•*•♦♦♦♦♦♦**•*♦♦******•♦♦*»*•••••**•••••*••»*•**•*• 

231  C  *  COMPUTE  THE  BEAM  WIDTH  BETWEEN  FIRST  NULLS.  * 

208  C  **♦•*••*♦•♦♦****♦♦*•* ♦♦♦I*********** ***♦♦♦*♦ ♦♦•»»**♦»** •****••••*•* 

208  REAL  U(NE).X(RE),T(>E).KII.KTI.K1P,KTP 

207  RAD= .174B3283E-01 

208  TP=,628318BE+01  -  -  -  . 

209  TM1K=0.0 

210  TMAI=12.0 

211  P  =  0.0 

212  NT- 140 

213  DT=(TMAX-TMIN)/NT 

211  C  »♦•*»*•*♦••*♦•»♦••♦♦♦♦•*«**••••*•**♦•♦♦****•••••**•*•••***•******** 

215  C  *  FIND  THE  MAXIMUM  VALUE.  * 

218  C  •*****♦**♦**••*•••*♦**•••*•*••••**•*•***»**•*♦***♦*•♦*••»**»**♦*••» 

217  CALL  PEKPOU(NE,0,I,T.KIP.KTP,PEAK) 

218  C 

219  C  *  RECOMPUTE  AND  NOW  FIND  THE  POINT  WHICH  IS  AT  THE  SLL  BELOW  THE  • 

220  C  *  PEAK.  • 

221  C  ***•♦♦•♦••****♦**•**♦•**♦***♦*♦••*••*•••*»*•••*•**♦•*•**•*••*•**•*• 

222  DO  3  I-0,NT 

223  T=TMIN+I»DT 

221  TPS=TP*SIK(RAD*T) 

225  KII=TPS 

228  CALL  P0WI(NE,U,X,T,KXI,0.0,P0W) 

227  IF  (POW/PEAK  .CE.  O.B)  THEN 

228  WRITE(2,100)  T,10.0*AL0G(P0W/PEAK) 

320  GO  TO  00 

230  ELSE 

231  ENU  IF 

232  3  CONTINUE 

233  100  F0R.MAT(F10.B,1X,F10,6) 

231  00  CONTINUE 

235  BWFN=2.0*T 

236  RETURN 

237  END 

238  C 

230  SUBROUTINE  KS?ACE( NMP ,NL ,DS ,SLL,HBAR,BBF ,KI ,KT ) 

210  C  •♦.••*•••**»•••***«•*•••••*•<*•••»••*•••*•»**•♦*****••♦♦»******.*** 

211  C  •  THIS  SUBROUTINE  SAMPLES  THE  KI  >0  REGION  OF  K  SPACE.  * 

212  C  **♦••**•*•*»***•******••*•*• •♦•••**♦•**•**•*♦••**•*•*♦**»********** 

213  REAL*-!  KX(*)  ,KT(*)  .KC,KR,KII,KTT,KXMIN,KIMAI,KTMIt:,KTMAI,KNULL 

211  RAD-.174B3203E-01 

215  PI- .3141B03E+01 


6G 


TP=.(J283186E+01 


24« 

247  C 

248  C 
24S  C 
280 
2E1 
262 

263 

264 
266 
266  ' 

267 

268 

269 
260 
261 
262 

263 

264  C 
266  C 
2«6  C 

287 
266 
260  C 

270  C 

271  C 

272  C 

273 

274  C 
276  C 

276  C 

277 

278 
2  78 
280 
281 
282 

283 

284 
28f 
386 

287 

288 
289 
29fi 

291 

292 

293 


*  COMPOTE  THE  FIRST  ZERO  OF  THE  THEORETICAL  BATLISS  PATTERS.  * 


S=ABS(SLL) 

C1=3.0387B3E-01 

C2=E.042622E-02 

C3=-2.798eE-04 

C4=3.43E-0e 

CB=-2.0E-8 

A=C1+S*(C2+S*(C3+S*(C4+S*CB))) 

C1=9.8B8302E-01 

C2=3.3388BE-02 

C3=l .40e4E-04 

C4=-l .9E-09 

C6=l .OE-08 

ZETA1-C1+S*(C2+S*(C3+S*(C4*S*C6))) 

01  =  (HBAR+0 . 6)  *SQRT(2ETAl*2ETAl/(  A*A4-HBAR)  ) 


•  COMPOTE  THE  K  VALOE  OF  THE  FIRST  BOLL. 


R=NL*DS 

KNULL=U1/R 


*  COMPOTE  THE  MAXIMOM  K  VALOE  ALLOWED  FOR  THE  FIRST  BOLL.  • 


KC=TP*0.4B 

KC=BBF*KNULL 


*  GENERATE  THE  RECTAROOLAR  LATTICE  OF  SAMPLE  POINTS. 


KTMIH=-TP 

KTMAX=TP 

KXMIN^O.O 

KXMAI=TP 

NKX=27 

NKT=46 

DKX=(KXHAX-KXHIN)/(NKX-1) 

DKY=(KT«AX-KTMIH)/(RKT-1) 

NMP  =  0 

DO  1  1=0, NKT 
KTT=KTMIN+I*DKT 
DO  2  J=O.BKX 
KXX=KXMIN+J*DKX 
KR=SORT(KIX*KXX+KTT*KTT) 

IF  ((KR  .CT.  KC)  .AND.  (KR  .LT.  TP) 

A  .AND.  (KXX  .GT.  0.0))  THEN 
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NHP=HMP+1 

KX(NMP)=K1X 
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2PS 

]B« 

397 

366  2 

IBS  1 

300 

301 
303  C 

303 

304  C 
SOB  C 
SOfl  C 

307  C 

308  C 
308 

310 

311  C 
313  C 

313  C 

314 

315 
318 

317 

318 
318 

320 

321 

322 

323 

324 

326 
328 

327 

328 
328 

330 

331 

332 

333 

334  C 

335  C 
338  C 

337 

338  C 

339  C 

340  C 
34  1 

342 

343 


KT(HMP)=KTT 
ELSE 
END  IF 
CONTINUE 
CONTINUE 
RETURN 
END 

SUBROUTINE  TAPER(KE.U.I.T.SLL.NBAR.BL,DS) 


*  INITIALIZE  THE  AMPLITUDE  DISTRIBUTION.  * 

*  NOTE  THAT  THE  AKPLITUDE=U«U.  THE  BATLISS  DISTRIBUTION  FORMULAS  * 

*  ARE  TAKEN  FROM  MODERN  ANTENNA  DESIGN  BT  T.  MILLIGAN.  * 


REAL*4  D(NE) ,I(HE) .T(HE) ,R00T(20) ,ZETA(4) ,DZER(20' ,fl(20) 
P1=.3141E93E+01 


*  DEFINE  THE  ROOTS. 


* 


R00T(l)=0.6860fl70 
ROOT(2)=1 .8C70608 
R00T(3)=2 .7171938 
ROOT(4)=3 .7281370 
R00T(6)=4. 7312271 
RG0T(6)=6. 7346206 
R00T(7):8.T38B281 
R00T(8)=7. 7386360 
R00T(9)=8. 7388606 
ROOT(10)=9. 740896 
R00T(11)=10. 7417436 
ROOT(12)=11 .7424476 
ROOT  I"  13)  =  12. 7430408 
root;  14)-.^13. 7436477 
ROOT(16)=14. 7439860 
ROOT(lO)=16. 7443079 
aOOJ (17)=16. 7447044 
ROOK  18)  17.7460030 
tOOT( 19}=ie .7462897 
ROOT(.  )--19. 7466093 


'  COMPOTE  APPROXIMATE  ARRAT  RADIUS. 


RADIU3=DS*NL 


•  CHOOSE  SIDELOBE  LEVEL. 


CALI.  CC1EFF(SLL,A,ZETA) 

CALL  NEVZ(A,ZETA,NBAR,ROOT,DZER) 
CALL  FBCOF(ROOT,OZER,NBAR,B) 


6^ 


344  C 

343  C 

343  C 

347 

348 

343 

330 

331 

33} 

333 

334 

*  COMPinE  AMPLITTOE. 

DO  1  K=1,NE 

R=SQRT(I(K)8X(K)+T(K)8T(K)) 

C0SPH1=I(K)/R 

SOM=0.0 

DO  2  M=1,HBAR 

RBO=R/RADIOS 

ARG=PI*ROOT(M)«RHO 

CALL  Jl( ARC, BESS 1) 

• 

3TB 

“  ~ S^H— 5UH^B\M} ^BESSl 

338  2 

CONTINUE 

367 

AMP=COSPHI*SOM 

3B8 

0(K)=AMP 

33R  1 

CONTINUE 

360 

CALL  NORHAL(BE,U) 

36] 

DO  3  1=1, NE 

38} 

WRITF,(2.*)  I. KD.KD.R. RADIOS. D(I) 

383  3 

CONTINUE 

384  100 

F0RMiT(IB,lI,F10.4.ir,F10.4,lX.F10.4) 

365 

RETURN 

see 

END 

387  C 

366 

SUBROUTINE  CO£FF(SLL.A,ZETA) 

399 

REiL44  ZETA<4> 

370 

S=ABS(SLL) 

371 

C1=3,0387G3E-01 

372 

C2=B.042e22E-02 

373 

C3=  2.r«e9E-04 

374 

C4S3.43E-00 

376 

C6=-2.0E-8 

376 

A=C14S»(C24S*(C3+S*(C4+S*CB))) 

377 

C1=B.8B8302E-01 

376 

C2=3.338e6E-02 

379 

C3=l .40fl4E-04 

360 

C4=-l .BE'OO 

361 

C6=l .OE-08 

362 

ZETA(1)=C1+S*(C2+S4(C34-S*(C4+S*C6)); 

363 

01=2.00337487 

384 

C2=1.141648E-02 

365 

C3=4.169E-04 

366 

C4=-3.73E-0« 

367 

C6=l .OE-08 

368 

2ETA(2)=Cl+S*(C2tS*(C3-S*{C4+S*C&))) 

360 

01=3.00638321 

390 

02=6.833945-03 

301 

03=2 .92eiE-04 

302 

C4=-l .81E-06 

69 


303 

394 

305 

390 

307 

396 

300 

400 

401  C 

402 

403 

404 

405 

406 

407  1 
406 
40P 

410 

411  2 
413 
413 

4M  C 

416 

416 

417 
416 
410 
■120 
•121 

422 

423 

424 

425 

426  2 

427 

428 
420 

4.10 

4.11 
432 
4.1.1 
4.14 

435  3 

4.16 

4.17 

436 
UP  1 
4  1(» 


ZET>(3)=C1+S*{C2+S*(C3+S*C4)) 

C1=4.00B18432 

C2=B.0179BE-03 

C3=2.173BE-04 

C4=-8.8E-07 

ZET1(4)=C1+S*(C24-S*(C3+S*C4)) 

RETURN 

END 

SUBROUTINE  NEUZ(A.Z£Ti,NBlR.aOOT»UZER) 
RE1L’*4  ZETA(4).R00T(20),UZER(20) 

“  DO  1  5=1,4  -  -  - 

ARG=ZETA(N)*ZETA(K)/(A*A'i-56AR*5BiR) 
OZER(H)=ROOT(HBAR+1)*SQRT(ARG) 
CONTINUE 
DO  2  5=B.NBAR-1 

ARG=(A*A+B*N)/(A*A+NBAR*RBAR) 

UZER(N)=R00T(MBAR+1)*SQRT(ARG) 

CONTINUE 

RETURN 

END 

SUBROUTINE  FBCOF(ROOT.UZER,HBAR,B) 
REAL*4  R00T(20),UZER(20),B(20) 

REAL*8  RAT,UN2,RM2,RN2,PR1,PR2 
Pl=3. 141693 
DO  1  M=1,NBAR 

RM2=R00T(M)*RC0T(M) 

PR1=1 .0 

DO  2  H=l,NBAR-l 

DN2=UZER(N)*UZER(N) 

RAT=1 .0-RM2/DN2 
PR1=PR1»RAT 
CONTINUE 
PR2-1 .0 
DO  3  N=1,NBAR 

IF  (N  .NE.  M)  THEN 
RN2=R00T(N)*R00T(N) 

RAT=1 .0-RM2/RN2 
PR2-PR2*RAT 
ELSE 
END  IF 
CONTINUE 
ARG=PI*ROOT(M) 

CALL  JKARG.BESSl) 
B(M)=(RM2/BESS1)*(PR1/PR2) 

CONTINUE 

RETURN 


4ti  END 


70 


442  C 

143  SUBROUTINE  J1(X,BESS1} 

444  C******************************* 

445  C  THIS  SUBROUTINE  COMPUTES  JIU) 

446  C******************************* 


c 

c 

'C 


447  10=-0. 60249086 

448  *1=0.21003673 

440  12=-0. 03964280 

ISO  *3=0.00443310 

451  *4=-0. 00031701 

452  *6=0.00001100 

453  80=0.70788460 

454  81=0.00000166 

455  82=0.01669007 

456  83=0.00017106 

457  B4=-0. 00240611 

458  86=0.00113063 

469  B6=-0. 00020033 

460  C0=-2. 36010440 

461  €1=0.12409612 

462  02=0.00006660 

46.3  C3  =  -0  00037879 

464  04=0.00074348 

465  CB=0. 00079824 

466  C0=-0. 00020160 

467  If  ((X  .GE.  0.0)  .*ND.  (X  .LE.  3.0))  THEN 

468  S=(X/3.0)*(X/3.0) 

469  P1=S*(A0+S*(*1+S*(*2+S*(*3+S*(*4+*B*S))))) 

470  BESS1=I6(0.6+P1) 

471  ELSE  IF  (I  .CT.  3.0)  THEN 

472  T=3.0/X 

473  Ql=B0+T»(Bl+T*(B2+T*(B3+T*(B44T*(BB+Be*T))))) 

474  Vl=ItC0+T»(Cl+T*(C2+T»(C34T*(C4+T*(CB+C6*T))))) 

475  BESS1=Q1*C0S(V1 )/SQRT(l) 

476  ELSE 

477  end  if 

476  RETURN 

4  79  END 

480  C 

481  SUBROUTINE  CGR*D(NE ,N5P ,NMP ,HS .HI .U.X ,T ,5 , HP , IP ,KI,KT ,GN ,H .P.Q ,R 

482  *,Z) 

483  C  ••*•»*»  ****************4**»***«4i*********»******»********4i**»*44*** 

481  C  ♦  COMPUTE  POWER  IN  *  SINGLE  DIRECTION.  * 

485  C 

486  REAL*4  U(NE) , I(NE) , T(HE) , S(RE) ,GB(HE) ,H(RE ,NE) , P (NE) ,Q(HE) ,R(HE) 

48'  REAL*4  Z(NE),KI(NSP),KT(NSP),KXI,KTI,KIP,KTP 

488  INTEGER  IP(HSP) 

489  C  ***+••«*•*♦•**♦•*•**•♦•**••*•*••**•**»****♦**•***************»«**** 

490  C  •  INITIALIZE  THE  NUMBER  OF  SIDELOBE  SAMPLE  POINTS.  • 


71 


4»l  C 

492 

493  C 

494  C 

495  C 
49e 

497  C 

498  C 

499  C 
soo  C 

501 

502  C 

503  C 

504  C 

505 

506 

507  C 

608  C 

509  C 

510  C 

511  C 

512 

613 

614 

515  2 

516  C 
51,-  C 
618  C 

519 

520 

521 

522 

SJ3  3 

524 

525 
528 

527  200 

528 

529  C 
630  C 

531  C 

532 

533  C 
531  C 

535  C 
539 
537 

536 
539 


HP=0 


♦  PERFORM  HS  RE-STIRTS.  * 


DO  i  L=1.NS 


*  FIND  THE  PEiK  SIDELOBE  POWER  POIHT  AID  UPDATE  TBS  HUMBER  * 

•  OF  MATCH  POUTS.  • 


CALL  PEAKSL(HE .O.I.T .HSP .HNP.KX.KT , IP .HP.PMAX) 


•  FIHD  THE  POWER  IH  THE  MAIH  BEAM.  • 


CALL  PEKP0W(HE,D,X,T,KXP,KTP,PE1K) 
PO^PEAK 


•  COMPUTE  AND  STOKE  AN  ARRAT  OF  THE  1ST  DERIVATIVE  OF  THE  ♦ 

•  POWER  IN  THE  MAIN  BEAM  DIRECTION  WITH  RESPECT  TO  THE  MTH  * 

•  variable.  * 


DO  2  M=1,NE 

CALL  DP0WI(HE.D.I.T.M.KXP.KTP,DPWI) 
S(M)=DPHI 
CONTINUE 


*  COMPUTE  AND  STORE  THE  AVERAGE  SIDELOBE  POWER. 


PA  =  0 .0 
DO  3  1=1. NP 

CALL  POWl(NE,U.X.T.KX(IP(I)),KT(IP(I)),POWERI) 
PA^PA+PUWERI 
CONTINUE 
PA=PA/NP 

WRITE(*,200)  L,10.0*ALOC10(PA/P0) ,10.0*ALOC10(PMAX/P0) 
WRITE(2,200)  L,10.0*ALOC10(PA/P0) ,10.0*AL0G10(PMAI/P0) 
F0RMAT(H,I6,1I,’AVE  PEAK  SLL  ’  .  II  ,F16 . 8  ,  II , 

*  'PEAK  SLL' ,11, FIB. 6) 


♦  COMPUTE  THE  NEGATIVE  OF  THE  GRADIENT. 


* 


DO  4  M=1,NE 


♦  COMPUTE  THE  DERIVA''IVE  OF  THE  AVERAGE  SIDELOBE  POWER.  • 


DPAM=0.0 
DO  B  1=1, HP 

CALL  DPOWI(NE.D,l,T,M,KX(IP(I)).KT(IP(r)) ,D?WI) 
DPAM=DPAM4DPWI 


72 


B40  6 
B4l 
S4l  C 

543  C 

544  C 
645 

B4e  4 

647  C 

648  C 

649  C 
660 
661 
652 
6J3 

654 

655 
568 
557 
658  8 

650 

sen 

sai 

662  7 

503  6 
664  C 
res  C 
566  C 
607  C 
668  C 
see  C 
h70 

571 

572 

573  12 

574  C 

575  C 

576  C 

577 

578 
570 
58(1 

561 

682  13 

583  C 

584  C 

585  C 
6Hfl 
687  C 
589  C 


COHTim 

DPiM=DPiM/VP 


*  COMPOTE  THE  lEGiTlTE  CRAOIEMT. 


* 


(}H(M)=(P1*S(M)/P0-DP1M)/P0 

CDRTIHUE 


•  COMPOTE  THE  HESSIAH.  * 


DO  0  M=1.KE 
DO  7  H=M,HE 

CALL  DDPOUKME.O.X.T.M.R.KIP.KTP.DDPWI) 

DDPO=DDPWI 
DDPA--0.0 
DO  8  1=1, NP 

CALL  DDPOW1(NE.O,I,7,M.»,KI(IP(I)),KT(IP(I)),DDPWI) 
DDPA=DDPA+DDPUI 
CONTIBOE 
DDPA=DDPA/HP 

H(M,N)=(GN(B)*S(M)+G»(M)*S(N)+DDPA-PA*DDPO/PO)/PO 

H(N,M)=H(M,H) 

CONTIBOE 

CONTIBOE 


•  START  THE  CONJUGATE  GRADIENT  ALGORITHM.  • 


*  initialize  the  residual. 


DO  12  M-1,RE 
P.(K)=GN(M) 
Z(M)=0.0 
COHTINOE 


•  INITIALIZE  SEARCH  VECTOR. 


* 


CALL  MATVEC(NE.H,R.Q) 
CALL  B0RH22(HE,Q,qN) 
BEOri .o/QN 
DO  13  1=1, BE 
P(I)=BEO*'Q(I) 
CONTINUE 


-  PERFORM  CONJUGATE  GRADIENT  ITERATIONS. 


DO  14  K=1,NI 


•  UPDATE  AMPLITUDE  VECTOR  AND  RESIDUAL. 


* 


73 


500 

c 

t*4i 

600 

CALL  M1TYBC(SE,H,P,Q) 

»et 

CALL  irORN22(HE,q.QN) 

602 

AK=1 .O/QH 

« 

603 

DO  IE  1  =  1. HE 

504 

Z(I)=Z(1)+AK*P(I) 

605 

R(I)=R(1)-AK*Q(I) 

506 

IB 

CONTINUE 

• 

507 

CALL  N0RM22(NE,R.RN) 

608 

ERR=SQRT(RN) 

500 

WR1TE(*.B00)  L.K.ERR 

eoo 

WRITE(2.B00)  L.K.ERR 

. 

601 

600 

FORMATOI. 'RESIDUAL' .  11.13. 11. 13. 11. ElB. 8) 

602 

C 

603 

C 

•  UPDATE  SEARCH  VECTOR. 

604 

C 

•  • 

60'. 

CALL  MATVEC(!;E,H.R.q) 

006 

CALL  N0RH22(NE,q,QN) 

007 

BEK=1 .0/QN 

608 

DO  16  1=1, HE 

000 

P(I)=P(I)+BEK«q(I) 

010 

16 

CONTINUE 

a:  i 

14 

CONTINUE 

612 

C 

•  •**•«**«•*******•.**«.•*•***•**••*.•**••**•*•*«,*»•*••*•*»**>.•» 

613 

C 

*  UPDATE  AMPLITUDE  VECTOR. 

t 

61  1 

C 

•  • 

015 

DO  IT  1  =  1, NE 

610 

U(I)=U(I)+Z(I) 

017 

17 

CONTINUE 

• 

6!6 

CALL  NORMAL(NE,U) 

010 

C 

***.***•«**«***•, •,««*,*,*,*****«*.,«**,«**4i,*«**,**«»*******«i 

620 

C 

*  WRITE  OUT  THE  NEW  ESTIMATE  TO  A  FILE.  * 

02! 

c 

•  • 

622 

REWIND  4 

623 

WRITE(4,101)  (I, 0(1). 1=1. NE) 

024 

101 

F0RKAT(IB,F1B.7) 

625 

IF  (NP  .EQ.  NMP)  THEN 

620 

GO  TO  99 

027 

ELSE 

628 

END  IF 

020 

1 

CONTINUE 

630 

00 

CONTINUE 

631 

RETURN 

0.^2 

END 

633 

c 

634 

SUBROUTINE  PEAKSL( NE ,U , I , T ,NSP ,NMP ,KX ,KT , IP , NP , PHAX ) 

630 

c 

•  THIS  SUBROUTINE  FINDS  THE  PEAK  SIDELOBE  POWER  LEVEL  COORDINATE 

• 

037 

c 

•  POII.’T  (KI(LKAX)  ,KT(LMAX)) .  IF  THE  POINT  IS  MEW  THEN  IT  IS 

* 

74 

# 


636  C 
630  C 
641)  C 

641 

642 

643  C 

644  C 

645  C 

646 
64  7 
648 
640 
060 
66 
053 
653 
664 

655 
650  1 
057  C 

656  C 
650  C 
oa  j  C 

601  C 

602  C 
663  C 
601 
d6C- 
060 
667 
606 
660 
070  2 
6V  I 
072 
673 
671 

0  7  5 
67('. 
677 
»178  C 
0/0 

r 

6Hl  C 

t'-H?  C 

6*.l 

OH'. 


*  ACCnMULATEI)  INTO  THE  POIHTER  ARRAT  IP.  THE  HHMBER  OF  POIHTS  IB  ♦ 

•  THE  IHTEGER  POIHTER  ARRAY  IP  IS  IP  SUCH  THAT  IP(I) ,1=1 . . .HP.  * 


REAL  0(NE) ,1(HE),T(HE)  KX{HSP)  ,KT(SSP) 
IHTEGER  IP(HSP) 


»  COMPOTE  THE  POWER  AT  EACH  POIHT  AHD  FIHD  THE  MU 


CALL  P0WI(HE,O,I.T,KX(l).KY(i),PWl) 

PMAX=PUI 

LMAX=1 

DO  1  1=1, HMP 

CALL  POWI(KE,0,X,T,KI(I),KT(I),PW1) 
IF  (pvr  GT.  PMAI)  THE.H 
LMAX=I 
PMAX=PWZ 
ELSE 
END  IF 
CONTINUE 


»  COMPARE  LK4X  TO  ALL  THE  PREVIOOS  POINTERS  IP(I)  ,1=1 . . .HP . 

»  IF  LMAI  MATCHES  ONE  OF  THE  PREVIOOS  POINTERS  THEN  THE 

•  POINTER  AW, AT  NEED  NOT  BE  UPDATED  (IFLAG  =  l).  IF  LMAX  DOES 

•  NOT  MATCH  A  PREVIOOS  POINTED.  VALUE  (IFLAG=0)  THEN  IT 

•  BECOMES  THE  HP+ITH  F-LEMEHT  OF  THE  POINTER  ARRAT . 


IFLAG=0 
2  1=1, NP 

IF  (IP(I)  .EQ.  LMAX)  THF.N 
IFLAC=) 

ELSE 
END  IF 
CONTINUE 

IF  (IFLAG  .Eq.  0)  THEN 
NP^^^NP*! 

IP(NP)=LMAX 
ELSE 
END  TF 
RETURN 
r.^-D 

SUBROUTINE  POWZ ( HE . D , 1 . Y .XX I . KTi , POWERI ) 


♦  COMPUTE  PO-'ER  IN  A  SINGLE  DIRECTION. 


REAL  U(NE) ,I(NE),y(NE) ,K1I,KTI 
SUM  y . 0 
no  1  K=1,!IE 
P.Sl^KXI^tCD'.TTI  +  KK) 


7ri 


987 

988  1 

989 
090 
991 
997  C 
903 
994  C 
099  C 
899  C 

997  C 

998 

999  C 

700  C 

701  C 

702 

703 
701 
705 
709  1 

707  C 

708  C 

709  C 

710 

711 

712 
'13 
7M  C 
715 
719  C 

717  C 

718  C 

719  C 

720 

721  C 

722  C 

723  t' 
721 

725 

726  C 

727  C 

728  C 
720 

730 

731 

732  C 

733 
731  C 
735  C 


S0M=SDW+0(K)*SIH(PSI) 

COKTINtfE 

POWERl=SaM*SOK 

RETORN 

END 

SUBROOTINE  DPOWKKE.O.I.T.H.KIl.KTl.DPWI) 


♦  COMPUTE  1ST  DERITITIVE  OF  POWER  WITH  RESPECT  TO  THE  MTH  VARIABLE* 

•  IN  THE  ITH  DIRECTION.  * 


REAL  0(NE),I(NE),T(NE).KII.KTI 


*  COMPOTE  TI  AND  VI. 


S0M  =  0.0 
DO  1  K-l.NE 
PSI=KII*I(K)+KTI*T(K) 
SaM=S0M+0(K)*SIN(PSI) 
CONTINUE 


♦  COMPUTE  1ST  DERIVATIVE  OF  TI  AMD  VI. 


PSIM=KII*I(M)+Kyi*T(M) 

DPWI=2.0*SIN(PSIM)*S0M 

RETURN 

END 

SUBROUTINE  DDPOWI ( NE ,0 , I , T ,M , N ,KXI ,KTI .DDPWI) 


•  COMPOTE  2ND  DEPH  ^  OF  POWER  WITH  RESPECT  TO  THE  M  AMD  HTH  • 

•  VARIABLES  IH  Tf  jIRECTION.  * 


REAL*4  U(NE) .7(«E) ,T(NE) ,KII ,KTI 


COMPUTE  2ND  DERIVATIVE  OF  TI  AMD  VI  WITH  RESPECT  TO  M,N.  * 


r5I«^KII*I(M)^KTI*T(M) 

PSIN=KII*I(N)*KTI*T(N) 


♦  COMPUTE  2ND  DERIVATIVE  OF  POWER. 


DDPWI=2. 0*SIH(P5IM)*SIH(FSIR) 

RETURN 

END 

SUBROUTINE  nOP.M22(  N  ,  A ,  AN) 


♦  THIS  SUBROUTINE  COMPOTES  THE  EUCLIDIAN  NORM  SQUARED  OF  A .  • 


76 


73«  C‘ 

737 

73a 

739 

740 

741  1 

742 

743 
74;'.  C 
7V' 

746  C 

747  C 
746  C 
740 

750 

751 

752 

753 

?h4  2 

735 

756  1 

757 
756 

7C0  C 

760 

76-  C 
7tJ2  C 
7(f*3  C 

7(^5 
766 
rev 
/ee 
;  60 
7/0 

771  1 

772 

773 

774  2 
TV5 

7  76 

777  C 

778 
770 

780  C 

781  C 
78?  C 
70.1 
784 


REAL*4  !(♦) 

AH=0.0 
DO  1  1  =  1, N 
AH=AN+A(I)*A(I) 

CONTINUE 

RETURN 

END 

SCBROUTXNE  MiTTEC(K,H.P,Q) 


•  THIS  SOBROOTIHB  COMPUTES  THE  MATRIX  PRODUCT  Q=>HP. 


REAL*4  H(r.N).P(K).q(N) 
DO  1  1^-1, H 
SUM^O.O 
DO  2  J=1 .N 
SUM  =  SUy.-H(I, .!)'■?(  J) 
rOKTINUE 
Q(I)=SU(1 
CONTINUE 
RETURN 
RIID 


SUBROUTIKE  MORMAL(H,A) 

•  *««»n**4t«*4.**«***««««»«*u**»*i********lt»*»»****«« 

•  TL'IE  SUBROUTINE  NORMALIZES  A  VECTOR  OF  LENGTH  B . 


REAL‘4  A(;;) 

VnAI  =  A(l ) 

DC  1  I-]  ,H 

xr  (A(l)  .CT.  VKAX)  TEEN 
VMAX^Ad) 

ELSE 
END  IF 
CONTINUE 
DO  2  i-i,:; 

A(I)  =  Aa)/ABS(VKAX) 
CONTINUE 
RETURN 
END 


SUBROUTINE  PEKPOV( NE ,U  . I , T  ,Kir ,KTP .PEAK) 


•  COMPUTE  POWER  IN  A  SIKJLE  DIRECTION  GIVEN  A  THETA  AND  PHI. 


REAL  U(NK) ,X(N£) ,T(NE) ,KII ,XTI .KIP.KTP 
RAD  -  .  174B,3203E-01 


/  t 


766 

TP=.«283186E*01 

T8« 

TKIH-0.0 

7sr 

TMAI=30.0 

766 

PMIK=-10.0 

760 

PMAX=10.0 

760 

sr=eE 

791 

»P=46 

703 

DT=(TM1I-TMIH)/*T 

703 

DP=(PMAX-  PMii;)/jrp 

794 

DO  1  1=0, 5T 

796 

T=T«IN+I*DT 

796 

TPS=TP*SIH(RAD*T) 

767 

DO  2  J=0,HP 

796 

P=PHIN+J*D? 

799 

RP=RAD*P 

600 

KXI=TPS'*COS(RP) 

601 

KT1=TPS*SIKCRP) 

802 

CALL  POUI(KE.O,I,T.KX1.KTI.POU) 

803 

IF  (I  .Eq.  1)  THEN 

804 

PEAK=POW 

«05 

Kir=KX7 

soe 

KTP^KTl 

807 

THEPK^T 

806 

PHI.nX-P 

806 

ELSE  IF  (POU  .CT.  PEAK)  THEN 

810 

PEAKt^POW 

81  1 

KIP--KII 

at: 

KTP-KTI 

813 

THEPK=T 

814 

PHIPK=P 

8U 

ELSE 

813 

END  IF 

8'7  2 

COKTIPUE 

«1A  1 

CONTINUE 

819 

WHITER • , * )  THEPK ,PHIPK ,KXP ,KTP .PEAK 

820 

WRITE(2,*)  THEPK, PBIPK , KIP ,KTr , PEAK 

821 

RETURN 

822 

END 

873  C 

624 

SUBROUTINE  CAIN (NE ,U ,I , T , CDB ) 

825  C 

8Jfl  C 

♦  COMPUTE  POWER  IN  A  SIKCLE  DIRECTION, 

• 

627  C 

828 

REAL  U(NF),I(nE),T(NE) ,KIP,KTP 

82(1 

RAD=.174B3293E-01 

830 

PI-.3141b0^E^Ol 

8.1 1 

TP=  .028318BEt01 

f32  C 

a.i:i  C 

•  GRAPH  INPUTS 

I* 

78 


83S 

ltP=41 

630 

NR=3S 

83? 

PMIM-0.0 

838 

pmi-300.0 

830 

TMIN=0.0 

840 

THAI=0O.O 

DP=CPMAI-PMIN)/(HP-1) 

842 

DR=(TMAI-T«IB)/(Sll  1) 

843  C 

9**19* 

€44  C 

*  COMPUTE  PEAK  VALUE. 

* 

84  S  C 

9**** 

840 

CALL  ?EKPOW(HE,0,I,T,KIP.KTP,PEAK) 

SI'  C 

9  *  <•  *  * 

846  C 

*  COMPUTE  THE  SOLID  ANGLE. 

* 

846  C 

950 

SUM^O.O 

8;.; 

DO  I  I^0,BR-2 

852 

T=7MIN+DR/2.0+I«DR 

853 

ST=SIN(RAD*T) 

8S4 

DO  2  J^0,»P-2 

855 

P=PMIN+DP/2.0+J*DP 

650 

CALL  power(i;e,u,i,t,t,p.pow) 

957 

SUM=SUH4ST*P0W 

ssa  2 

cohtihoe 

8611  1 

COHTItlOE 

HHI) 

SOLID  ^RAD*RAD*DR4-DP*S0M/PEAK 

fifll 

CDB^10,0*AL0G10(TP/S0LID) 

902 

retur:: 

803 

EN'D 

894  C 

805 

SUBROUTl!i£  POWPATCNE  ,U ,  I  .T) 

8  19  C 

»**** 

8ilV  C 

•*  cntli'OTF.  POWER  Iri  A  SIWGLE  DIRECTION . 

* 

r.Ps  C 

»**** 

800 

HEAL  U(.V£)  ,I(}IE)  ,T(KE)  .KIP.KTP 

P70 

OPf;N<  L'NIi -  7  ,  F J  LE^-  ’  DATA  ’ ) 

971 

RAD=.  174&3293i'  01 

872 

PI-  .3r41L93E  +  0.'. 

67:5 

TP--  .028318BE401 

874  C 

|i**** 

I|75  C 

8/0  C 

•  GRAPH  IKPUTS 

* 

977 

KP-ill 

9  70 

NR=  20 

879 

NS-O 

980 

PMIN=0 .0 

89| 

PMAX=3eo  0 

887 

ThlN^O .0 

79 


683 

TMAI=B0.0 

684 

DP  =  (PMAI-PmH)/(MP-l) 

68B 

DR=(TMiI-TMIF)/(ini-l) 

686 

AX=1B0.0 

687 

AT=AI 

886 

AZ^flSO.O 

889 

PO=6S.O 

890 

qo=6o.o 

69t 

o 

II 

803 

IB=2.0 

803 

TA=-2.0 

804 

TB=2.0 

695 

ZA=-2.0 

606 

ZB=2.3 

607 

AT=.3 

8P8 

WRITE(7,*)  HP,HR,HS.AI,AT,AZ 

800 

WRITE(7,.)  iA.ifl,Tl.TB,2A,ZB 

000 

WRITE(7,*)  AT.PO.qO 

HOI  C 

902  C 

•  PHI  VALUES. 

903  C 

004 

DO  1  I=0,NP-1 

005 

P=PHIN4I*DP 

006 

WRITE(7.*)  RAD*P 

907  1 

COHTINOE 

909  C 

909  C 

*  RADIAL  (THETA)  VALUES. 

* 

910  C 

01  1 

DO  2  I=0,NR-1 

912 

T=TM1N+I*DR 

913 

WRITE(7,*)  BAD*T 

9)4  2 

COHTINUE 

91S  C 

Pit)  C 

•  FUNCTION  VALUES. 

917  C 

916  C 

919  C 

*  COMPUTE  PEAK  VALUE. 

* 

920  C 

031 

CALL  PEKPOW(HE,U,J,T,KIP,KTP,PEAK) 

922  C 

923  C 

♦  WRITE  OUT  THE  FORMALIZED  POWER. 

921  C 

025 

DO  6  1=0, NR- 1 

020 

T=TMIN+I*DR 

927 

DO  8  J=0,NP-1 

026 

P=PMIN+J*DP 

029 

CALL  POWER(NE,U,I.T,T,P,POW) 

030 

WRITE(7,*)  POW/PEAK 

931  6 

CONTINUE 

80 


031  B 

CONTIRVE 

933 

RETTOH 

931 

EKD 

035  C 

036 

SOBROUTIHE  SPOWPAT(«E,0,1,T) 

037  C 

033  C 

*  COMPUTE  POWER  IR  i  SIRGLE  DIRECTION. 

• 

030  C 

040 

REAL  U(HE),X(RE).T(RE) 

041 

OPEN  (ONIT-^e  ,F1LE= '  SDATA  • ) 

043 

RAD=.174B3293B-01 

043 

PI-.314IB93E+01_ _ 

044 

TP=.0283186E+01 

B45  C 

046  C 

*  GRAPH  INPUTS 

9 

1147  C 

99*999 

v/4e 

HP=119 

040 

NR-B3 

050 

NS-0 

9&1 

PMIN=0 .0 

083 

PMAX=3flO.O 

0E3 

TMIN=0.0 

084 

THAX=90.0 

085 

DP=(PMAX-PMIN)/(NP-1) 

066 

DR=(TMAX-TMIN)/(NR-1) 

067 

AX=flOO.O 

088 

AT=AX 

080 

AZ=AI 

eeo 

PO-6b.O 

901 

qOrflO.O 

062 

iA=  2.0 

003 

IB=2.0 

064 

T4--2 .0 

048 

TB--2 ,0 

008 

7,A---2.0 

967 

ZB=2.3 

008 

AT^.3 

0(30 

h:*'i 

WIIITE(8,0  RV, NR, NS, II. AX. AT, A^ 

071 

WRITE(8,')  XA.IB.TA.TB.ZA.ZB 

072 

WRITEO,*)  AT.PO.QO 

i»73  C 

0  74  C 

*  PHI  VALUES. 

0 

075  C 

9  78 

DO  1  I=0,KP-1 

0  77 

P-P«INH*DP 

0  78 

WRITE(8,*:i  R1D*P 

970  i 

CONTINUE 

esn  C 

81 


961 

c 

*  RiOlAL  (THBTi)  TALtJBS. 

• 

SSI 

c 

963 

DO  2  I=0,IR-1 

•64 

T=TMIB+I*DR 

986 

WRITE(8.*)  IUD*T 

960 

2 

COHTIBtrB 

967 

C 

SB9 

c 

•  niBCTIOH  VALDES. 

980 

c 

SBO 

c 

CSl 

c 

*  COMmS  PBAX  VALOE. 

903 

c 

993 

PBAK=- 100.0 

604 

DO  3  I-O.BR-l 

S65 

T=TM1JI+I*DR 

ooe 

DO  4  J=0,HP-1 

697 

P=PM1H+J»DP 

9ea 

CALL  P0WER(HE,O,I,T,T.P.P0W) 

606 

IF  (POW  ,CT.  PEAK)  THEH 

1000 

PEAX=POW 

1001 

ELSE 

1002 

ESD  IF 

1003 

4 

COSTINDE 

1004 

3 

COKTINDE 

1005 

C 

looa 

C 

*  WRITE  ODT  TBE  NORMALIZED  POWER. 

• 

1007 

C 

*»«** 

1008 

DO  6  I=0,NR-1 

looo 

T=TMIH+I*DR 

1010 

DO  e  J=0,*P-1 

101 1 

P=PMIB+J*DP 

1012 

CALL  P0WER(1E,D,X,T,T,P,P0W) 

1013 

WRITE(8,*)  POW/PEAK 

1014 

e 

CONTINOE 

1016 

6 

CONTIBDE 

loie 

RETORN 

1017 

END 

loie 

c 

1019 

SDBRODTINE  POWER( SE ,D ,I , T ,T ,P ,POU) 

\070 

c 

1021 

c 

*  COMPUTE  POWER  IN  X  SINGLE  DIRECTION  GIVEN  A  TBETA  AND  PBI. 

0 

1022 

c 

10?  t 

REAL  U(NE),X(NE),T'NE),KXI,KTI 

1024 

RAD:=  .ir4B32B3E-01 

l'>26 

TP=.e28318BE+0l 

lU2n 

TPS=TP*SIN(RAD*T) 

1027 

RP=RAD*P 

102* 

KXI=TPS*CDS(RP) 

1029 

KTI=TPS*SIN(RP) 

82 


1031 

DO  1 

1032 

PSI=KX1*I(K)+KTI*T(K) 

1033 

S0M=S0M+D(K)'>S1H(PSI) 

1034  1 

COHTIB0E 

t03S 

POW=SnM*SOM 

1036 

RETURK 

1037 

END 

